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fatigue  crack  growth  rata  (In. /cycle). 

any  discrete  value  of  fatigue  crack  growth  rate  (in. /cycle), 
maximum  fatigue  crack  growth  rate  (in. /cycle), 
minimum  fatigue  crack  growth  rate  (in. /cycle), 
inverse  fatigue  crack  growth  rate  (cycles/in.). 

Napierian  base  (2.718281828). 
e raised  to  the  x powar. 

expected  frequencies  in  the  chi>square  test. 

cumulative  density  function  of  corrected  data. 

inverse  cumulative  density  function  for  the  generalized  4>para* 
meter  ganma  distribution. 

density  function  of  a random  variable. 

standard  normal  probability  function. 

shape/power  parameter  for  the  gamma  distributions. 

Golden  Section  search  constant  (0.618033989). 

ganma  probability  fun  cion. 
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vairable  used  in  ueriving  C . 
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variable  used  in  deriving  C . 

iteration  counter  for  the  interior  point  penalty  function 
algorithm. 
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variable  used  in  deriving  C . 

number  of  equiprobable  intervals  for  the  chi-square  teat. 

maximum  likelihood  estimator  function. 

natural  logarithm  (base  a). 

logarithm  (base  10). 

slope. 

cumulative  load  cycle  count. 
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final  cuonilatlve  load  cycle  count. 

any  discrete  cumulative  load  cycle  coimt. 

average  cumulative  load  cycle  count  used  by  the  secant  method. 

log  (base  10)  scaled  cycle  count  data  used  by  the  log* log 
Incremental  polynomial  methods. 

scaled  cycle  count  data  used  by  the  Incremental  polynomial 
methods. 

number  of  data  points  In  a data  set. 

number  of  data  points  lost  at  each  end  of  the  data  by  the 
incremental  polynomial  methods. 

number  of  distribution  parame te¥f^^^ 

number  of  data  points  in  Che  incremar^ed  strip  in  the 
incremental  polynomial  methods. 


observed  frequencies  in  the  chi-square  test.\^ 

objective  function  used  by  the  interior  point  penalty  function 
algorithm. 

maximum  applied  load  during  the  load  cycle  (lbs.), 
minimum  applied  load  during  the  load  cycle  (lbs.), 
chi-square  test  class  end  points. 

R ratio. 

coefficient  of  multiple  determination. 

iteration  variable  in  the  interior  point  penalty  function 
algorithm. 

standard  deviation  of  Che  errors. 


standard  deviation  of  the  standard  deviations. 


standard  error. 


SSRES  residual  sum  of  squares. 

T location  parameter  in  the  Golden  Section  search  method. 

TCSS  total  corrected  sum  of  squares. 
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variable  In  the  gamma  probability  function. 

variable  used  in  deriving  da/dN  for  the  quadratic  log-log 
7-polnt  Incremental  polynomial  method. 

variable  of  Integration  In  the  chi-square  tall  area  equation. 

covariance  matrix. 

inverse  covariance  matrix. 

variable  to  be  plotted  on  the  abscissa. 

mean  of  the  errors. 

mean  of  the  standard  deviations. 

variable  to  be  plotted  on  the  ordinate. 

Kolmogorov- Smirnov  teat  statistic, 
standarized  distribution  variate, 
power  parameter  In  the  gamma  distributions, 
hypotheals  acceptance  level. 

shape  parameter  in  the  log  normal  distributions, 
gamma  function. 

Euler's  constant  (0.5772157). 

distance  between  two  consecutive  half  crack  length  data  points 
(In.  or  mm. ) . 

charge  in  stress  Intensity  during  the  load  cycle  (ksi-square 
root  In.). 

any  discrete  value  of  AK  (ksl-square  root  In.). 

number  of  cumulative  load  cycles  between  two  consecutive 
cumulative  load  cycle  count  data  points. 

change  In  cumulative  load  cycle  count  normalized  by  the  change 
In  half  crack  length  (cycles/ln.) . 

change  In  applied  load  during  the  load  cycle  (lbs.). 

convergence  criterion  constant. 

positive  constant  approaching  zero. 
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e characteristic  value  of  the  3-paraineter  Welbull  distribution.  j 

scale  parameter  (sometimes  the  mean)  of  a distribution. 
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V degrees  of  freedom. 

TT  PI  (3.141592654). 

a shape  parameter  (sometimes  the  standard  deviation)  of  a 

distribution. 

T location  parameter  (sometimes  the  terminus)  of  a distribution.  \ 

X random  variable. 
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X chi-square  test  statistic. 

X^  corrected  random  variable  data. 

X^  any  discrete  value  of  a random  variable. 
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SECTION  1 


INTRODUCTION 

Throughout  the  course  of  history.  It  has  always  been  desirable  to 
be  able  to  predict  the  life  of  a given  design  under  expected  service  con- 
ditions. Life  prediction  in  natal  structures  has  necessitated  a need  for 
knowledge  about  the  sKtal  fatigue  phenomenon.  The  natal  fatigue  process, 
as  it  is  known  today,  is  conplex  and  is  still  not  fully  understood.  There 
are  nany  variables  which  influence  the  life  of  a metal  structure,  such  as 
the  sMterial,  loading,  and  geometric  characteristics  of  the  particular 
structure.  This  investigation  involves  only  the  determination  of  the  ef* 
feet  of  Qiaterial  properties  on  life  prediction. 

One  of  the  primary  mechanisms  by  which  metal  fatigue  occurs  is  the 
propagation  of  microscopic  cracks  [ij.  The  study  of  fatigue  crack  propa- 
gation behavior  has  been  widely  conducted  for  some  time  in  an  effort  to 
understand  metal  fatigue  more  fully.  The  information  obtained  from  crack 
propagation  studies  is  then  used  in  estimating  the  fatigue  life  of  struc- 
tures and  components.  Ideally,  it  is  desirable  that  this  estimated  life 
will  exactly  predict  the  actual  life.  Unfortunately,  there  are  many  vari- 
ables which  Influence  this  prediction  and  some  are  not  well  understood. 

Ons  of  the  most  important  of  these  variables  is  how  well  the  empirical 
crack  growth  relationships  obtained  from  experimental  data  actually  re- 
present the  observed  crack  propagation  behavior. 

The  raw  data  from  a fatigue  crack  propagation  test  are  the  half 
crack  length,  a,  and  the  number  of  cumulative  load  cycles,  N,  needed  to 


grow  Che  crack  Co  ChaC  lengch  from  some  reference  InlClal  crack  lengch 
for  sllghcly  Increasing  stress  incenslCy  level  load  conditions,  called 
constant  amplitude  loading.  A plot  of  typical  raw  fatigue  crack  propa- 
gation data  is  shown  in  Figure  1.  The  current  interpretation  of  this 
raw  data  focuses  upon  the  fatigue  crack  growth  rate  as  a function  of  an 
applied  stress  intensity  parameter,  usually  AK,  the  change  of  the  stress 
intensity  during  the  load  cycle.  The  fatigue  crack  growth  race  Is  de- 
fined as  the  rate  of  extension  of  Che  crack  with  respect  to  the  number 
of  applied  load  cycles  [2].  Actual  deteminatlon  of  the  crack  growth 
rate  requires  an  evaluation  of  the  slope  of  the  raw  a vs.  N data  at 
various  discrete  points,  which  results  in  Che  derivative  of  a with  re- 
spect to  N,  normally  called  da/dN.  A plot  of  typical  da/dN  vs.  AK  data 
is  shown  in  Figure  2. 

The  importance  of  Che  fatigue  crack  growth  rate  as  a variable  of 
interest  is  bom  out  in  the  fact  that  the  fatigue  crack  growth  rate  is 
nearly  independent  of  the  geometry  for  the  same  stress  intensity  level 
of  loading  [3].  This  allows  crack  growth  behavior  prediction  based  only 


on  the  knowledge  of  the  crack  growth  rate  vs.  the  stress  intensity  level  ; 

of  loading  for  a given  material  for  any  geometry  chosen.  Obviously,  this 

I 

! 

would  be  an  important  design  tool  if  the  crack  growth  behavior  predictions 
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were  accurate  and  reliable.  These  crack  growth  behavior  predictions  can 
be  used  to  predict  the  number  of  load  cycles  needed  to  grow  a crack  from 
an  initial  crack  length,  a^,  to  soiae  new  crack  length,  a^,  and  the  dis- 
tance a crack  propagates,  Aa,  during  a specified  number  of  applied  load 
cycles.  In  addition,  using  various  prediction  techniques,  the  constant 
amplitude  loading  crack  growth  rate  behavior  is  used  to  predict  variable 
amplitude  loading  crack  growth  rate  behavior  [4] . 
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There  are  several  nethods  of  numerically  determining  the  crack 
growth  rate  from  the  raw  a va.  N data.  It  has  been  suspected  that  the 
crack  growth  rate  calculation  method  has  a very  significant  effect  on 
the  variance  of  the  resulting  growth  rate  va.  stress  intensity  parameter 
data  [2,5,6,^ . 

During  the  prediction  of  crack  growth  behavior,  the  crack  growth 
rate  vs.  stress  intensity  parameter  data  is  integrated  back  to  obtain 
predicted  a vs.  N behavior.  Considerable  variation  in  this  predicted 
crack  growth  behavior  has  been  experienced,  thus  hindering  accurate  life 
estimates  [2, 5, 6, 7].  This  variation  is  a result  of  variation  in  the  raw 
crack  growth  data,  variation  due  to  the  crack  growth  rate  calculation 
method,  and  material  variations. 

This  investigation  will  compare  several  numerical  growth  rate  cal- 
culation methods  and  attempt  to  find  the  method  which  introduces  the 
least  amount  of  error  into  the  growth  rate  vs.  stress  intensity  parameter 
data.  It  will  also  attempt  to  describe  crack  growth  behavior  in  a sta- 
tistical manner  with  the  expectation  that  this  statistical  description 
of  crack  growth  behavior  will  reduce  the  large  amount  of  error  currently 
present  in  life  prediction. 


5 


SECTI(»J  11 


BACKGROUND 

Mecal  fatigue  has  long  been  recognized  as  a random  phenomenon  [8j, 
but  until  recently,  little  effort  was  devoted  to  applying  statistical 
tools  to  fatigue  crack  propagation  behavior.  By  fitting  different  equa- 
tions to  the  crack  growth  rate  vs.  stress  intensity  parameter  data,  nu- 
merous equations  of  fatigue  crack  growth  have  been  suggested  [9].  How- 
ever, due  to  scatter  In  the  data,  it  has  been  impossible  to  select  which 
equation  Is  the  most  appropriate.  Also,  when  the  original  crack  growth 
data  are  predicted  from  these  equations,  the  correlation  with  the  original 
data  Is  generally  very  poor  [81.  Due  to  the  large  amount  of  scatter  in 
the  crack  growth  rate  vs.  stress  intensity  parameter  data,  investigators 
have  started  using  statistical  methods  to  characterize  fatigue  crack  pro- 
pagation behavior  [5, 6, 7, 8].  It  can  be  easily  shown  that  the  amount  of 
data  scatter  is  generally  considerably  greater  than  can  be  accounted  for 
by  experimental  inaccuracies  [8j.  It  has  been  pointed  out  that  the  re- 
maining scatter  is  due  to  the  essentially  random  nature  of  fatigue  crack 
growth  which  is  a result  of  the  relative  nonhomogrnelty  of  the  material 
[8,10]. 

From  a macroscopic  viewpoint,  it  is  often  convenient  to  regard  a 
metallic  material  as  a homogeneous  continuum,  and  basing  engineering 
calculations  on  this  assumption  does  not  generally  lead  to  serious  error. 
However,  the  scatter  observed  in  fatigue  testing  of  a metallic  material 
arises  precisely  because  it  is  not  a homogeneous  continuum,  when 
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considered  on  a microscopic  scale  [8],  Consequently,  it  is  important  to 
examine  fatigue  crack  growth  from  a statistical  viewpoint.  In  order  to 
include  fatigue  crack  propagation  scatter  in  the  general  overall  charac* 
terisation  of  fatigue  crack  propagation  behavior,  this  investigation 
will  apply  statistical  concepts  to  fatigue  crack  growth  behavior. 

In  considering  the  crack  growth  from  some  Initial  crack  length,  a^, 
to  a new  crack  length,  a^,  there  is  a certain  mean  and  variance  associated 
with  the  number  of  load  cycles  required  for  this  amount  of  crack  growth 
which  characterizes  the  statistical  distribution  of  N at  a^.  A schematic 
representation  of  this  distribution  of  N is  shown  in  Figure  3.  In  order 
to  statistically  characterize  the  crack  growth  behavior,  it  is  necessary 
to  determine  the  distribution  of  N from  experimental  tests. 

The  variance  In  N Illustrated  above  can  be  due  to  random  errors  in 
the  measurement  of  a,  N,  and  AK,  to  systematic  errors  in  these  measure- 
ments, and  to  the  statistical  variation  in  the  material's  growth  rate 
properties.  Through  the  use  of  accurate  equipment,  the  random  errors  in 
the  measurement  of  a,  N,  and  AK  cst  be  reduced  to  an  acceptable  level  and 
measured  by  a separate  test.  Through  a careful  experimental  set  up  and 
procedure,  the  systematic  measurement  error  can  be  reduced.  From  this, 
the  desired  statistical  behavior  of  the  material's  crack  growth  proper- 
ties can  be  determined. 

In  considering  the  crack  growth  rate  vs.  stress  intensity  parameter 
data,  there  is  some  statistical  distribution  associated  with  the  crack 
growth  rate,  da/dN,  at  some  stress  intensity  level,  AK^.  A schematic 
representation  of  this  distribution  of  da/dN  is  shown  in  Figure  4.  In 
order  to  statistically  characterize  the  crack  growth  rate  behavior,  it 
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is  also  necessary  to  determine  the  distribution  of  the  crack  growth  rate 
from  experimental  tests. 

The  variance  of  da/dN  illustrated  above  originates  in  the  variance 
present  in  the  original  a vs.  N data.  The  density  of  the  raw  data  (es> 
sentially,  the  distance  between  2 consecutive  data  points,  &a)  and  the 
crack  growth  rate  calculation  method  both  contribute  to  the  overall  vari- 
ance of  da/dN.  In  order  to  determine  the  variance  of  da/dN  due  to  the 
variance  in  the  original  a vs.  N data,  it  is  necessary  to  determine  the 
effect  of  both  data  density  and  the  crack  growth  rate  calculation  method 
on  the  variance  of  da/dN. 

Once  the  crack  growth  rate  vs.  stress  intensity  parameter  data  has 
been  obtained,  the  next  step  is  to  be  able  to  predict  the  change  in 
crack  length  for  a given  number  of  applied  load  cycles  or,  Inversely,  the 
number  of  applied  load  cycles  for  a given  change  In  crack  length.  The 
variance  of  this  prediction  is  directly  related  to  the  variance  of  the 
crack  growth  race.  In  order  to  evaluate  the  effectiveness  of  this  pre- 
diction, it  is  necessary  to  predict  the  original  a vs.  N data  from  the 
crack  growth  rate  data  and  then  compare  the  predicted  a vs.  N data  with 
the  original  a ve.  N data. 

This  a vs.  N prediction  can  be  accomplished  by  either  of  two  methods. 
The  currently  popular  method  is  to  numerically  integrate  the  mean  da/dN 
vs.  dK  curve  to  obtain  predicted  a ve.  N data  [2, 4, 9, 6, 7, 9] . Howevar, 
no  adequate  method  for  determining  the  resulting  scatter  in  a or  N axiste 
[5].  An  alternate  method  uses  the  knowledge  of  the  distribution  of  da/dN 
and  the  fact  that  da/dN  is  an  independent  random  variable  to  obtain  a ve. 

N step  by  step.  This  method  is  discussed  in  detail  in  Section  7.3.  Using 
this  method,  the  variances  of  both  a and  N can  be  readily  obtained. 
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SECTION  III 


OBJECTIVES  OF  IMVESTIOAIION 


The  cuin  purpose  of  this  investigation  was  to  apply  statistical  con* 
cepts  and  theory  to  the  study  of  fatigue  crack  propagation  behavior.  In 
doing  this,  there  were  four  main  objectives  to  be  net.  They  were; 

1)  Determine  the  statistical  distribution  of  N (cuDulatlvo 
load  cycle  count)  as  a function  of  a (crack  length). 

2)  Determine  which  crack  growth  rate  calculation  method 
yields  the  least  amount  of  error  when  the  crack  growth 
rate  curve  is  integrated  back  to  the  original  a vs.  H 
data. 

3)  Determine  the  statistical  distribution  of  da/dN  (crack 
growth  rate)  as  a function  of  AK  (stress  intensity 
parens  ter. 

A)  Determine  the  variance  of  a set  of  a vs.  N data  predlcceo 
from  the  da/dN  distribution  parameters. 
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SBCTIOM  IV 


CRACK  GROWTH  RATH  CALCUUTION  METHODS 

Nunaroua  aathodf  of  calculating  Cha  crack  growth  rata  from  tha  raw 
a va.  N data  have  been  ueed  by  varioui  inveatigatore  [2,3,6].  None  of 
theae  aeem  to  be  unlvereally  accepted,  but  rather  each  Invaatlgator 
aeemad  to  favor  a different  nethod.  Since  it  wee  virtually  lapoesible  to 
inveatigete  ell  of  theae  nethode,  eix  of  the  oiore  important  methoda  were 
aeleeted  for  examination.  Thaee  methoda  are: 

1)  The  aecent  method, 

2)  The  modified  eecent  method, 

3)  The  linear  7*point  incremental  polynomial  method, 

4)  The  quedretic  7-point  incremental  polynomial  method, 

5)  The  linear  log- log  7*polnt  Incremontal  polynomial  method,  end 

6)  The  quadratic  log-log  7-point  incremental  polynomial  method. 

4.1  Secant  Method 

The  eecent  method  ie  a finite  difference  nethod  and  perhepe  the 
eimpleet  of  the  nethode  coneidered  [2,3,6] . Baeicellyr  the  accent 
method  celculatee  the  elope  of  a etreight  line  between  2 adjacent  a ve. 

N data  pointe.  It  then  epproximatee  thie  elope  ee  the  elope  of  the 
tangent  line  of  the  a ve.  N curve  at  an  average  crack  length,  e^,  end 
average  cycle  count,  M^.  A echematic  reproeentetion  of  thn  eecent 
method  le  ehown  in  Figure  3. 
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(CRACK  LENGTH) 


The  average  crack  length,  a^,  la  given  by 


r - 

*1  2 


Similarly,  the  average  cycle  count,  la  given  by 


IT  . 


(1) 


(2) 


The  slope  of  the  line  connecting  the  2 adjacent  data  points,  which  Is 
used  to  approximate  the  growth  rate,  la  given  by 


da 

dN, 


<‘fH  • ‘l> 


at  a^  and 


(3) 


4.2  Hod if led  Secant  Method 

The  modified  secant  method  la  really  an  extension  of  the  secant 
method.  Basically,  this  method  averages  the  growth  rates  obtained  by 
Che  secant  method  so  that  Che  da/dN  data  coincides  with  the  original  a 
va.  N data.  The  beginning  and  end  points  are  assumed  to  be  equal  to  the 
first  and  last  growth  rates,  respectively.  A schematic  representation  of 
Che  modified  secant  method  is  shown  In  Figure  6. 

The  growth  rate  is  given  by 


(4> 


at  a^  and  for  1«2  to  (n-1)  where  n la  the  number  of  data  points  in  the 
data  set. 

The  first  growth  rate  data  point  la  given  by 


at  a^  and 


da 

dN 


^*2  • 


V 


1 


(N,  - Nj) 


(5) 
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n (CROCK  LENGTH! 


The  laet  growth  rate  data  polot  ie  given  by 


da  *D-1^ 


at  a and  N . 
n n 


4.3  Linear  7»Polnt  Incremental  Polynomial  Method 
The  linear  7-polnt  Increnental  polynomial  method  la  the  simplest  of 
the  four  incremental  polynomial  methods.  In  each  of  the  incremental 
polynomial  methods,  a polynomial  is  fit  by  the  method  of  least  squares 
to  a series  of  data  points t called  a strip,  and  the  derivative  of  the 
polynomial  is  evaluated  at  the  middle  point  [2,5,6].  This  strip  is 
then  incremented  by  one  data  point  and  the  curve  fitting  and  evaluation 
process  is  repeated.  The  strip  incrementation  process  is  repeated  until 
all  of  the  data  points  have  been  used.  Any  odd  number  of  data  points 
can  be  used  for  the  incremented  strip,  although  7 points  are  usually 
used.  The  incremental  polynomial  leethods  differ  basically  in  the  poly* 
nomial  which  is  fit  to  the  data. 

Initially  the  strip  data  points  are  scaled  in  the  following  manner. 
Two  constants,  and  C^,  are  calculated  as  follows  [5,6]; 


N,  - N, 

^~"h8 


where 


n , - I 

strip 


Where  n . is  the  number  of  data  points  in  the  strip.  Note  that  C.  is 
strip  * 

the  center  of  the  strip  cycle  count  data  and  is  the  range  of  the  strip 
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cycle  count  data.  The  data  scaling  is  then  performed  as  follows: 


N 


S 


\ : 

(10) 


whore  is  the  scaled  cycle  count  data.  As  a result  of  this  scaling, 
the  strip  cycle  count  data  runs  from  -1  to  +1.  This  insures  that  when 
least  squares  curve  fitting  occurs,  the  scale  of  the  data  will  not  in- 
fluence the  curve  fitting,  which  is  a constant  danger  when  using  least 
squares  as  a curve  fitting  technique. 

After  the  curve  fitting  has  been  performed,  the  derivative  of  the 
resulting  polynomial  is  then  evaluated  at  the  midpoint  of  the  strip,  N^. 
This  evaluation  takes  into  account  Che  scaling  that  was  performed  prior 
to  the  curve  fitting.  A schematic  representation  of  the  incremental 
polynomial  method  is  shown  In  Figure  7. 

In  Che  linear  7-point  incremental  poljrnomial  method,  the  fitted 
polynoisial  is  s first  order  linear  straight  line.  After  fitting  by 
linear  least  squares,  the  fitted  polynomial  takes  the  following  form: 


Substituting  the  scaling  equation. 


(11) 


(12) 


Taking  the  derivative  of  a with  respect  to  N, 


(13) 


Obviously,  for  a straight  line,  the  slope  is  independent  of  where  the 
derivative  is  evaluated  at. 
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R (CRRCK  LENGTH) 


n (CRRCK  LENGTH) 


4, A Quadratic  7-Point  Incremental  Polynomial  Method 
The  quadratic  7*polnt  Incremental  polynomial  method  has  gained  wide 
acceptance  as  a valid  crack  growth  rate  calculation  oaethod  [5»6].  In 
this  method,  the  fitted  polynomial  is  a second  order  curve.  After  fitting 
by  second  order  least  squares,  the  fitted  polynomial  takes  the  following 
form; 


Substituting  the  scaling  equation, 


(W) 


(15) 


Taking  the  derivative  of  a with  respect  to  N and  evaluating  at  the  mid' 
point,  N^, 


(16) 


4,5  Linear  Log»Log  7-Polnt  Incremental  Polynomial  Method 

The  linear  log- log  7-polnt  Incremental  polynomial  method  was  used 
to  determine  If  Che  deCn  could  be  llnearleed  by  a log^^  transformation 
on  both  the  crack  length  and  cycle  count  data.  This  method  is  essentially 
the  same  as  the  linear  locrsmencal  polynomial  method  except  for  Che  log 
transformations  of  the  input  data  Just  prior  to  the  data  scaling. 

The  fitted  polynomial  Is  linear  and  takes  the  following  form 


lo,  . - b,  * 

Where  Is  the  leg  scaled  cycle  count  data. 


(17) 


4, A Quadratic  7-Polnt  Incremental  Polynoinlal  Method 
The  quadratic  y- point  incremental  polynomial  method  haa  gained  wide 
acceptance  as  a valid  crack  growth  rate  calculation  method  [5,6].  In 
this  method,  the  fitted  polynomial  Is  a second  order  curve.  After  fitting 
by  second  order  least  squares,  the  fitted  polynomial  takes  the  following 
form: 


a - b + b,N„  + b,N-‘ 
o IS  2 S 


Substituting  the  scaling  equation, 

2 


a ■ 


(14) 


Taking  the  derivative  of  a with  respect  to  N and  evaluating  at  the  mid- 
point, N^, 


(16) 


4.5  Linear  Log-Log  7-Polnt  Incremental  Polynomial  Method 

The  linear  log- log  7-point  Incremental  polynomial  method  was  used 
to  determine  if  the  data  could  be  linearized  by  a IoSj^q  transformation 
on  both  the  crack  length  and  cycle  count  data.  This  method  is  easentlally 
the  same  as  the  linear  Incremental  polynomial  method  except  for  the  log 
transformations  of  the  Input  data  Just  prior  to  the  data  scaling. 

The  fitted  polynomial  Is  linear  and  takes  the  following  form: 


log  a - b^  + b^Hj^ 

where  is  the  leg  seeled  cycle  count  data. 


(17) 
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The  derivation  of  this  equation  is  shown  in  Appendix  A. 


4.6  Quadratic  Loa-Loa  7-Point  Increaental  Polynomial  Method 

The  quadratic  log- log  7*point  incremental  polynoiaial  method  was  used 
to  determine  if  a second  order  curve  fit  could  improve  the  performance  of 
the  linear  log-log  7-poinc  incremental  polynomial  method.  This  method  is 
essentially  the  same  as  the  linear  log-log  7-point  incremental  polynomial 
method  except  that  the  fitted  polynomial  is  second  order  instead  of  first 


order. 


The  fitted  polynomial  takes  the  following  form: 


1.8  . - b,  * . bjH^' 


The  growth  rate,  da/dN,  for  this  method,  evaluated  at  the  midpoint,  N^, 


is  given  by 


b^  ^ _ 

L c,  ” c,  J . L 


2b^logN^  - 2b^C^ 


(logN^)^  - 2b^C^logN^  ^ 

10  ^ c? 


The  derivation  of  this  equation  Is  shown  in  Appendix  B. 
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SECTION  V 


STATISTICAL  CONCEPTS 


When  used  properly,  statistics  Is  extremely  useful  in  quantifying  the 
results  of  many  engineering  experiments.  In  many  applications,  however, 
statistics  is  used  as  a quick  substitute  for  a thorough  experimental 
analysis  and  often  tliaes  it  is  used  without  checking  the  underlying 
assumptions  or  else  the  results  are  misinterpreted.  In  an  attempt  to 
alleviate  these  problems,  the  statistical  concepts  used  in  this  invest!* 
gation  and  their  use  as  tools  in  analyzing  fatigue  crack  growth  behavior 
will  be  presented  and  discussed. 

5.1  Histograms 

The  first  step  in  statistically  analyzing  any  set  of  data  is  to  see 
what  the  data  looks  like.  Histograms  are  statistically  derived  pictures 
of  a data  set.  They  give  a rough  idea  of  the  shape  of  the  density  func* 
tion  of  the  data.  They  also  ^ive  a rough  estimate  of  the  average  value 
and  the  amount  of  variability  present  in  the  data. 

The  most  conson  histogram  used  is  a frequency  histogram.  The  data  is 
divided  into  several  classes  and  the  frequency  of  the  data  in  each  class 
is  plotted  against  the  limits  of  the  classes  [11].  This  type  of  histogram 
freqiiently  takes  the  form  of  a bar  chart.  A slight  modification  of  this 
involves  calculating  the  relative  frequencies  in  each  class  by  dividing 
the  frequency  in  each  class  by  the  total  number  of  data  points.  The  re- 
lative frequencies  are  then  plotted  against  the  limits  of  the  classes. 
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This  is  called  a relative  frequency  histogram  [12].  An  example  of  a re- 
latl\e  frequency  histogram  Is  shown  in  Figure  8. 

Another  convenient  form  of  the  histogram  Is  called  a cumulative  fre- 
quency histogram.  This  histogram  shows  the  frequency  of  data  less  than  or 
equal  to  a specified  value.  It  la  calculated  by  cumulatively  adding  suc- 
cessive class  frequencies  of  the  frequency  histogram  from  the  smallest 
class  value  to  the  largest  class  value.  It  frequently  takes  the  form  of 
a step  chart.  Again,  the  relative  cumulative  frequencies  can  be  calculated 
by  dividing  the  cumulative  frequencies  by  the  total  number  of  data  points 
so  that  the  last  value  of  the  relative  cumulative  frequency  is  equal  to 
one.  When  the  relative  cumulative  frequencies  are  plotted  against  the 
limits  of  the  classes,  the  resulting  plot  is  called  a relative  cumulative 
frequency  histogram  [12].  An  example  of  a relative  cumulative  frequency 
histogram  is  shown  in  Figure  9. 

5.2  Distributions 

Once  a rough  idea  of  what  the  density  function  of  the  data  looks  like 
based  on  the  histograms,  the  next  step  is  to  try  to  fit  the  data  to  several 
likely  distributions.  Eight  different  distributions  were  selected  as 
likely  candidates  for  the  distribution  of  fatigue  crack  propagation 
variables. 


S.2.a  Two-Parameter  Normal  Distribution 
The  most  widely  used  distribution  in  statistics  is  the  two-paranecer 
normal  distribution  [12].  This  distribution  was  selected  as  a candidate 
for  the  distribution  of  fatigue  crack  propagation  variables  mainly  for 
this  reason  and  for  the  sake  of  completeness. 
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Figure  8.  Typlcel  Relative  Frequency  Hiatogram 
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Figura  9.  Typical  Ralatlva  Cumulatlva  Fraquancy  Hlatograo 


The  two  pareneters  of  the  two-paraneter  normal  diatributlon  are  Che 
mean,  dealgnated  by  which  la  the  acalc  parameter,  and  the  atandard 
deviation,  daaignated  by  o,  which  la  Che  ahape  parameter.  The  denalty 
function,  f(x),  for  the  tvo'< parameter  norauil  diatributlon  la  given  by 

C 11,12]. 

1 r 1 - « < X < ■ 


“'{•2  ^ 
u ^ ^ ^ 


The  aatioiatea  for  the  mean  and  atandard  deviation  ace  computed  by  [11,12] 

n 

. _ iSi^ 


.ti  <’<1  • »' 


Where  n ia  the  number  of  data  poiota  and  the  aymbol  aymbollees  an 
eacimated  value. 

The  atandard  errora  of  Che  eacioatea  provide  a meaaure  of  how  good 
theae  eacimatea  are.  The  atandard  errora  of  the  eatlmated  aiean  and 
atandard  deviation  are  given  by  [13] 


S.  E.  a -/f 
I.  f.  5 • ^ [(n«l)  - 2.  / 


where  F repreaenta  the  gaama  function.  The  covariance  of  (I  with  S la 
alwaya  equal  to  sero,  due  to  their  orthogonality  [13]. 


5.2.b  Two'ParaMCer  Log  Normal  Diatributlon 
The  two'paraaneer  log  normal  diatributlon  haa  bean  auapected  of  being 
a likely  candidate  for  the  diatributlon  of  fatigue  crack  propagation 
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variables  [5, 6].  Essentially,  Che  two* parameter  log  normal  distribution 
states  that  the  log^Q  o£  Che  random  variable  x>  i>a>f  log^QXt  £•  normally 
distributed. 

The  two  parameters  of  the  cwo«paramecer  log  normal  distribution  are 
Che  scale  paranmter,  and  g,  Che  shape  parameter.  Tho  density  function 
for  Che  two'parometer  log  normal  distribution  Is  given  by  [14,15] 

h ^ 5 5 • 

o < p V • 

The  estimates  for  u *nd  g arc  computed  by  using  the  following  aqua* 
tions  [14] 


(27) 


e • i^i  - u.y 


(28) 


The  standard  errors  of  these  estimaCes  are  given  by  the  following 
equations  [13,14]. 


S 


s. 

. i 


(29) 


(30) 


The  covariance  of  ^ with  g Is  again  always  equal  to  sero,  due  to  their 
orthogonality  [13]. 


5.2.C  Three-Parameter  Log  Normal  Distribution 
With  the  expectation  of  a better  fit  of  Che  data,  the  three -parameter 
log  noroul  distribution  was  considered  as  a candidate  for  the  distribution 
of  fatigue  crack  propagation  variables.  The  main  difference  between 
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the  two- parameter  and  the  three-parameter  log  normal  distributions  Is  the 
Inclusion  of  the  location  parameter  In  Che  three -parameter  log  normal  dis- 
tribution. 

The  three  paraowtars  of  the  three -parameter  log  normal  dlatrlbtitlon 
are  the  scale  parameter,  g,  Che  shape  parameter,  and  the  terminus,  r, 
which  Is  the  location  parameter.  The  density  function  for  Che  three-para- 
meter log  normal  distribution  Is  given  by  [14,15] 


£(x) 


o-t')  /sy 


- «e  < < ce 

0 < g < OB 

- OB  < T < X 


(31) 


The  difficulty  In  using  distributions  containing  a location  parameter 
Is  the  estimation  of  that  location  parameter.  The  parameter  estimation 
methods  used  to  obtain  the  value  of  the  location  parameter  are  presented 
In  Section  5.3. 

Once  Che  location  paramscar,  t,  has  been  astlmated,  u and  g are  astl- 
maced  using  the  following  equations  [15]. 

n 


- T) 


(32) 


n 


To 

values, 
don  Is 


(33) 


n 

obtain  the  standard  errors  ot  the  estimates  and  the  covariance 
the  covariance  matrix  for  the  three-parameter  log  normal  dlstrlbu- 
computed.  The  covariance  matrix  is  a synmetric  matrix  and  Is  given 


by  [16] 
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(34) 


r 


I 

I 

I 


I 


r - - 

(i  + 1)  exp(p)  “ 2p 


V - C, 


where 


8 

“ 2 3 


.axp[?  - g 
2^(^  + 1)  ®xp(p)  - ij  2 p expj^;!  - g 

«xp[2ii  - g 


Cy  - -J* (35) 

n^(P  + 1)  exp(p)  - 2 p - ij 

The  etenderd  errors  of  the  estimates  are  given  by  the  diagonal  terns  and 
the  covariances  between  the  estimates  are  given  by  the  o£f>dlagonal  terms 
of  the  3 by  3 covariance  matrix. 


5.2.d  Three •Parameter  Welbull  Distribution 
The  three •parameter  Welbull  distribution  has  long  been  considered  In 
representing  fatigue  data  [17].  For  this  reason,  the  three^parameter 
Welbull  distribution  was  considered  as  a candidate  for  the  distribution  of 
fatigue  crack  propagation  variables.  This  distribution  also  Includes  the 
location  parameter  as  one  of  Its  three  parameters  and  thus  the  difficulty 
of  Its  estimation  arises.  Two  basic  methods  were  used  to  estimate  the 
parameters  (Section  5.3)  and  each  loethod  required  different  parameters. 
Thus,  two  secs  of  Welbull  parameters  and  their  associated  equations  will  be 
presented. 

The  first  set  of  the  three  parameters  of  the  three>paraneter  Welbull 
distribution  Include  the  characteristic  value,  9,  '.dilch  Is  the  scale  para^ 
meter,  the  Welbull  slope,  fi,  which  is  the  shape  paranatsr,  and  the  expected 
minimum  value  of  x>  which  Is  the  location  parameter.  The  density  func- 
tion for  these  parameters  Is  given  by  [18] 
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i 

i 
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In  Che  method  of  estiraaClng  the  location  parameter  used  with  this  set  of 
parameters,  all  three  parameters  are  estimated  simultaneously. 

The  second  set  of  three  parameters  of  the  three-parameter  Weibull  dis- 
tribution include  b,  the  scale  parameter,  c,  the  shape  parameter,  and  the 
terminus,  t,  which  is  the  location  parameter.  The  two  sets  of  parameters 
are  related  as  follows. 


b - 6 - X., 


c - B 


The  density  function  for  Che  second  sec  of  parameters  is  given  by  [19] 


f(x)  - cCx-Tl'^'^b'^ex  Ami  » 


0 < b < * 


“•<T<X 

As  with  the  previous  set  of  parameters,  all  three  parameters  are 
estimated  simultaneously  when  the  location  parameter  is  estimated.  To 
obtain  Che  standard  errors  of  the  estimates  obtained  by  the  method  re- 
ferred to  above  and  Che  covariance  values,  the  covariance  matrix  for  the 
three-parameter  Weibull  distribution  is  computed.  The  covariance  matrix 
is  a symmetric  matrix  and  is  given  by  [13] 


V - v‘ 


where 
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*2 

c 

f2 


£i£dLi  r(i-i) 

b c 


^ r(i4) 

b 


(42) 


where  y i®  Eulers  Constant  (0.577215)  and  ^ represents  the  diganma  func- 
tion. 

Once  the  covariance  matrix  is  obtained,  the  standard  errors  of  the 
estimates  and  the  covariances  between  the  estimates  are  obtained  from  the 
same  terms  in  the  covariance  matrix  as  outlined  above  for  the  three- para- 
meter log  normal  distribution. 


5. 2. a Ganna  Distribution 

Due  to  the  nature  of  the  fatigue  crack  propagation  process,  two  im- 
portant assumptions  can  be  made.  The  first  assumption,  called  the  in- 
creasing failure  rate  assumption,  states  that  b"  <tu8e  the  crack  growth 
rate  Increases  as  the  crack  grows  (under  consr  plitude  conditions), 

the  rate,  or  probability,  of  failure  increases  we  crack  grows.  The 
second  assumption  states  that  the  distribution  of  a fatigue  crack  propa- 
gation variable  is  independent  of  the  crack  length  and  is  a function  of 
the  Initial  crack  length  only.  If  these  two  assumptions  are  mede,  then 
it  can  be  proven  that  a generalised  gamma  distribution  is  a valid  distri 
bution  for  any  fatigue  crack  propagation  variable  [13,20], 
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Geperolited  Four- Parameter  Ganma  Distribution.  The  four  parameters 
of  the  generalized  four- parameter  gamma  distribution  are  the  location  para- 
meter, T,  the  power  parameter,  or,  the  scale  parameter,  b,  and  the  shape/ 
power  parameter,  g.  The  shape  parameter,  c,  is  given  simply  by  [21] 

c - g cr  (43) 

The  density  function  for  the  generalized  four-parameter  gamma  distribution 
Is  given  by  [21] 


b*  r(g)  bao 


(44) 


b 2 0 
g i 1 


All  four  parameters  are  estimated  simultaneously  using  the  parameter 
estimation  methods  presented  in  Section  5.3.  To  obtain  the  standard  er- 
ros  of  the  estisutes  and  the  covariances  between  the  estimates  the  covari- 
ance matrix  for  the  generalized  four- parameter  gamns  distribution  is  com- 
puted. The  covariance  matrix  is  a symnetrlc  matrix  and  Is  given  by 
[13,21.22] 


„ 1 -1 
V > — V 
n 


(45) 


where 


M— 


♦ '(I) 


l<  ^ g tiil] 

b 

tiAl 

a 


r(|4) 

a a 

T2 


b r(B) 


o V rd; 


r(|)’ 

[u'  -7i  . i] 

?r(i)  '■ 


(46) 


where  ^ ' represents  the  trlgamoia  function. 
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The  standard  errors  of  the  estimates  are  given  by  the  diagonal  terms 
and  the  covariances  between  the  estimates  are  given  by  the  off-diagonal 
terms  of  the  A by  4 covariance  matrix. 

Three-Parameter  Gamma  Distribution.  If  the  power  parameter,  a.  Is 
assumed  to  be  equal  to  one,  the  generalized  four-parameter  gamma  distri- 
bution reduces  to  Che  three-parameter  gamma  distribution.  The  density 
function  for  the  three-parameter  gamma  distribution  is  given  by  {21,23] 

/v_  r X>T 


b®  r (g)  *-  b J g > 1 


The  three  parameters  are  estimated  using  the  same  method  used  for  the 
generalized  four-parameter  gamma  distribution.  The  standard  errors  and 
covariances  are  found  by  using  the  covariance  matrix  for  the  generalized 
four-parameter  gamma  distribution  (equations  AS  and  46)  and  setting  a 
equal  to  one. 

Generalized  Three-Parameter  Gamma  Distribution.  If  the  fatigue  crack 
propagation  variable  of  interest  Is  ^NMa,  then  considering  the  fatigue 
crack  propagation  process  it  would  be  expected  that  AN  would  be  zero  for 
Aa  zero  {24J.  From  this,  It  Is  assumed  that  the  location  parameter,  y* 

Is  equal  to  zero  which  reduces  the  generalized  four-parameter  gamna 
distribution  to  the  generalized  three-parameter  gamma  distribution.  The 
density  function  for  the  generalized  three-parameter  gamma  distribution 


Is  thus  I 2 3] 


X > 0 


f(X)  » 


fa  bM  • b > 0 

r(g)  l-  ■> 
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The  standard  errors  of  the  estimates  are  given  by  the  diagonal  terns 
and  the  covariances  between  the  estimates  are  given  by  the  off-diagonal 
terms  of  the  4 by  4 covariance  matrix. 

Three-Parameter  Gamma  Distribution.  If  the  power  parameter,  a,  is 
assumed  to  be  equal  to  one,  the  generalized  four-parameter  gamma  distri- 
bution reduces  to  the  three-parameter  gamma  distribution.  The  density 
function  for  the  three-parameter  gamma  distribution  is  given  by  [21,23] 

r . . Ill  („) 

b*  r (g)  ^ b j 8 1^ 

The  three  parameters  are  estimated  using  the  same  method  used  for  the 
generalized  four-parameter  gamma  distribution.  The  standard  errors  and 
covariances  are  found  by  using  the  covariance  matrix  for  the  generalized 
four-parameter  gamma  distribution  (equations  45  and  46)  and  setting  a 
equal  to  one. 

Generalized  Three-Parameter  Gamma  Distribution.  If  the  fatigue  crack 
propagation  variable  of  interest  is  ^M/as,  then  considering  the  fatigue 
crack  propagation  process  it  would  be  expected  that  aN  would  be  zero  for 
Aa  zero  [24] . From  this,  it  is  assumed  that  the  location  parameter,  y, 
is  equal  to  zero  which  reduces  the  generalized  four-parameter  gatmna 
distribution  to  the  generalized  three-parameter  gaona  distribution.  The 
density  function  for  the  generalized  three-parameter  gamma  distribution 
is  thus  [23] 


f(X) 


b*“  r(s) 


X ^ 0 
0^1 
b > 0 
g > 1 


(49) 
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The  three  parameters  are  estimated  using  the  same  method  used  for  tin 


genaralited  four- parameter  gamma  distribution.  The  standard  errors  and 
covariances  are  found  by  using  the  3 by  3 submatrix  for  b,  g,  and  a from 
the  4 by  4 covariance  matrix  for  the  generalized  four-parameter  gamma  dis- 
tribution (equations  45  and  46). 

Two-Parameter  Ganna  Distribution.  If  the  power  parameter,  o,  is  again 
assumed  to  be  equal  to  one,  the  generalized  three -parameter  ganna  distribu- 
tion reduces  to  the  two-parameter  gaona  distribution.  The  density  function 
for  the  two-parameter  gamma  distribution  is  given  by  [11,  12,  23] 


f(x) 


. 1. 


g-1 


b«  r(g) 


exp[- 


X i 0 

, b 2 0 
• i 1 


(50) 


The  two  parameters  are  estimated  using  the  same  method  used  for  the  gener- 
alised four -parameter  gaotta  distribution.  The  standard  errors  and  covari- 
ances are  found  by  using  the  3 by  3 submatrix  used  for  the  generalized 
three -Parana ter  gamna  distribution  and  settings  equal  to  one. 


5.3  Parameter  Estimation  Methods 

Since  the  determination  of  the  estimates  of  the  parameters  is  critical 
to  a proper  fitting  of  the  data  to  the  two,  three,  and  four-parameter  dis- 
tributions, two  different  parameter  estimation  methods  were  used  [14,  23]. 
The  first  nethod,  a graphical  method,  wes  selected  for  its  simplicity 
[17,  18,  26,  27].  The  second  method,  the  method  of  maximum  likelihood 
estimators  (ML£),  was  selected  because  of  its  reliability,  accuracy,  and 
widespread  acceptance  [14,  15,  19,  21-23,  28-323. 

5. 3. a Graphical  Method 

This  method  was  the  first  of  the  two  methods  attempted,  due  mainly  to 
its  simplicity  in  use  [17,  18].  This  method  was  tried  with  both  the 
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three -parameter  log  normal  distribution  and  the  three-parameter  Welbull 
distribution.  The  graphical  method  Involves  plotting  the  data  on  special 
probability  paper  whose  axis  scales  correspond  to  special  distribution 
characteristics  and  then  selecting  the  value  of  the  location  parameter 
such  that  the  resulting  plot  of  data  follows  a straight  line  [17,181. 

Once  the  estimate  of  the  location  paratseter  is  known,  the  estimates  of 
the  other  two  parameters  are  made  graphically.  Since  only  three  para- 
meters can  be  estimated  graphically,  this  limits  the  use  of  this  method 
to  two  or  three  parameter  distributions  [27]. 

For  the  three-parameter  log  normal  distribution,  a plot  of  Y vs.  X 
yields  a straight  line  for  data  that  follows  a three-parameter  log  normal 
distribution  [14]  where 


Y “ G(z) 


(51) 


X - log^o(x^) 


(52) 


where  G(z)  is  the  equation  for  the  standard  normal  probability  scale  which 
is  given  by  [11] 

2. 


(53) 


where 


F(Xc) 


(54) 


where  F(x^)  is  the  cumulative  density  function  of  the  corrected  data, 
is  the  value  of  the  data  corrected  for  the  value  of  the  location  parameter 
by  the  following  equation. 


Xc  “ X - Xp 


(55) 
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For  the  three -parameter  Weibull  distribution,  a plot  of  Y vs.  X where 


(ryfe)) 


X - In  (x^)  (57) 

yields  a straight  line  for  data  that  follows  the  three-parameter  Weibull 
distribution  [17,183. 

For  both  of  these  plots,  F(x^)  corresponds  to  the  median  ranks  which 
are  calculated  by  [18] 


^ ^ ^c  " ^58) 

To  determine  the  value  of  the  location  parameter  such  that  the  re- 
sulting plot  yields  a straight  line,  an  iterative  process  which  minimizes 
some  variable  must  be  used.  For  the  graphical  method,  the  variable  to  be 
minimized  is  the  curvature  of  a second  order  curve  fit  using  least  squares, 
thereby  assuring  a straight  line.  One  of  the  fastest  and  most  efficient 
of  the  many  minimization  methods  available  is  the  Golden  Section  search 
method  [26] . 

In  the  Golden  Section  search  method,  the  value  of  the  curvature  (the 
variable  to  be  minimized)  is  calculated  at  two  optimal  locations  and, 
based  on  these  values,  a certain  area  where  the  curvature  minimum  is  knowri 
not  to  exist  is  excluded  from  the  rest  of  the  search.  This  process  is 
repeated  until  the  area  remaining  to  be  searched  is  less  than  some  toler- 
ance level.  The  value  of  the  location  parameter  in  this  area  is  then 
taken  as  the  estimated  value  of  the  location  parameter.  A schematic  re- 
presentation of  the  Golden  Section  search  method  is  shown  in  Figure  10. 
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( CURVfiTURE  ] 


T (LOCflTION  PflRPIMETER) 

IF  BgCXi) < BgCXg) 

THEN  Xg  BECOMES  X^ 

IF  BgtXj)  > BaCXg) 

THEN  Xi  BECOMES  X„in 

IF  Be(Xi)  = BglXg) 

THEN  Xe  BECOMES  X««( 
AND  Xi  BECOMES  X„in 

H r 0.618033989 

Figure  10.  Golden  Section  Seerch  Method 
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5.3.b  Maxloun  Likelihood  Betiiuttora  MeChod 
After  the  graphical  nethod  was  perfected  and  uc'^.d,  the  need  for  a 
more  statistical  approach  to  the  estimation  of  the  ^ raoeters  of  the  two, 
three,  and  four-parameter  distributions  became  evident  (Section  8.1). 

This  led  to  the  use  of  the  Maximum  Likelihood  Estimators  method  to  statis- 
tically estimate  the  distribution  parameters. 

The  Maximum  Likelihood  Estimators  (KLE)  method  involves  solving  maxi- 
mum likelihood  equations  through  the  use  of  a nonlinear  programming  algor- 
ithm [14,15,19,21,22,23,28,30,31,32].  Many  forms  of  the  maxioum  likelihood 
equations  have  been  determined  by  investigators  for  the  three -parameter  log 
normal  distribution,  the  three-parameter  Welbull  distribution,  and  the  two, 
three,  and  four- parameter  gamma  distributions  [ 14, 15, 19,26'.32]  . 

Three -Parameter  Log  Kormal  Distribution.  The  maximum  likelihood  equa- 
tion used  in  this  investigation  for  the  three- parameter  log  normal  distri- 
bution is  [15] 


In  L(t)  - - n[a(T)  + I In  8 (t)]  (59) 

tSi  tn(x^-T) 

where  ii(T)  ■ (60) 

n 

i-i  [^"  - ii(T)] 

and  8 (T)  • * (61) 

Three -Parameter  Welbull  Pistribuclon.  The  maximum  likelihood  aqua- 
tion used  in  this  investigation  for  the  three-parameter  Welbull  distribu- 
tion is  [19] 

L(b,c,T)  ■ n(ln  c - clnb)  (c  - 1)  E In(x.-T)  - b’^  E (62) 

i-l  ‘ 1-1  ^ 

Note  chat  the  maximum  likelihood  equation  Is  a function  of  all  three  para- 
omters  whereas  for  Che  three-parameter  log  normal  distribution,  the  maxioum 
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likelihood  equation  is  a function  of  Just  the  location  parameter.  However, 
the  scale  parameter,  b,  ran  be  factored  out  of  this  equation  and  estlTiated 
separately.  The  resulting.  < parameter  maximum  likelihood  equation  for 
the  three -parameter  WeibuH  distribution  is  [13] 


L(c,t)  - In  c - In  [i  In(Xi-T)] 


(63) 


where  the  estimation  of  the  scale  parameter  is  given  by  [13] 

j.. 


n 


,)'t 


(64) 


The  effect  of  reducing  the  number  of  parameters  in  the  maximum  likelihood 
equation  is  to  reduce  the  computing  time,  and  thus  the  cost,  of  the  maxi- 
mization of  the  maximum  likelihood  equation. 

Generalited  Four -Parameter  Gamma  Distribution.  The  maxionim  likelihood 
equation  for  Che  generalized  four-parameter  gamma  distribution  is  [21,30, 
32] 

n 
Z 

1-1 

(65) 

Inr(g) 


r n _ 

L(b,g/i',T)  - n Ino  + (go-l)  I I In(x.-T)  I - gn  ln(b®) 

^ J _1  ^ J 

- E [-i— J - n 
i-1  b 


The  number  of  parameters  in  this  equation  can  also  be  reduced  by  factoring 
out  the  scale  parameter,  b.  The  resulting  three  parameter  maximum  likeli- 
hood equation  for  the  generalized  four- parameter  gamma  distribution  is  [13] 

b<g,cr,T)  - n lno+  (gcr-l)  Z ln(Xj^-T)J  - gnj^l  + In  ^ Z (X^-t)®^  - IngnJ 

(66) 

- n Inr(g) 

where  the  estimation  of  the  scale  parameter  is  given  by  [13] 
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(67) 


.O'  l/or 


Three -Parameter  Gomna  Dlgtrlbution.  The  maximum  likelihood  equation 
for  the  three-parameter  gamma  dlatributlon  reduced  to  eliminate  the  scale 
parameter  is  [13] 

L(g,T)  - (g-1)  ^ ln(x^-T)J  - gn  ^1  + In  ^ 


(68) 


- n Inr(g) 

where  the  eatimatlon  of  the  scale  parameter  is  given  by  [131 


(69) 


Generalised  Three -Parana ter  Gamma  Distribution.  The  maximum  likeli- 
hood equation  for  the  generalised  three -parameter  gamma  distribution  Is 
[21,30,32] 


[n  _ 

E ln(x.)|  - gn  ln(b®) 
1-1  ^ 

- r [— T ‘ n Inr(g) 

1-1  ^ b-^ 


(70) 


This  equation  can  also  be  reduced  to  eliminate  Che  scale  parameter,  b. 
The  resulting  two  parameter  maximum  likelihood  equation  for  Che  gener- 
alized three -parameter  gamma  distribution  Is  [13] 

L(g/if)  - n linr  + (gor-l)  In(x^)]  - gn^l  ^^1^**^  ' ^"^8n)J 


(71) 


- n Inr(g) 

where  the  estimation  of  the  scale  parameter  Is  given  by  [13] 
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(72) 


Twp-Parmeter  Gaaiaa  Distribution.  The  maxlmun  likelihood  equation 
for  the  two'paraneter  gamna  distribution  reduced  to  eliminate  the  scale 
parameter  is  [13] 


L(g)  " (g-1)  [ 2^  ln(x^)  ]-■{'  + In  (X^)^  - ln(gn)J 

- n lnr{g) 

where  the  estimation  of  the  scale  parawter  is  given  by  [13J 


(73) 


(74) 


Interior  Point  Penalty  Function.  After  some  experience  using  the 
Graphical  method  to  estimate  the  location  parameter  of  the  three«para> 
meter  Weibull  distribution,  it  was  found  chat  Che  iteration  tended  Co  go 
to  minus  infinity  in  some  cases.  Since  this  was  the  global  (overall)  maxi- 
mum of  the  function  to  be  maximized,  it  became  necessary  to  use  a method 
that  converged  on  the  local  maximum,  and  not  the  global  maximum.  The 
method  used  to  achieve  this  requires  the  use  of  an  interior  point  penalty 
function  which  prevents  the  value  of  each  of  the  parameters  from  reaching 
either  of  its  global  limits  [IS]. 

The  interior  point  penalty  function,  better  known  as  the  objective 
function  in  nonlinear  programming  terms,  for  the  three -parameter  log  normal 
distribution  using  the  aiaximuffl  likelihood  equation  is  [IS] 


P(t 


,r)  - In  L(t)  - r^(T  + (x^n  ' t (75) 


25. 


where  c^  is  a large  positive  number  10  ), 
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r la  an  Iteration  variable,  and 

-8 

c_^  la  a email  positive  number  <^10  ). 

The  objective  function  for  the  three •parameter  Welbull  dlatrlbution 
ualng  the  two  parameter  maximum  likelihood  equation  la  [IS] 

P(T,c,r)  - L(t,c)  - ^(t  + "*■  ^\,tn  " 


■1 


(76) 


+ (C  - €^)  ‘ + (10  - C - 

The  objective  function  for  the  generalleed  four^paraneter  gaana  dla- 
trlbution using  the  three  parameter  maximum  likelihood  equation  Is  [IS] 

P(gia.T,r)  - L(g/»,T)  - r[(g  - 1 - a^)‘^  + (100  - g - 


+ to  - 1 - + (100  - a - t^)'^  (77) 

♦ <T  + ♦ <X^„  - T - 

The  objective  function  for  the  three -parameter  gamma  distribution 
using  the  two  parameter  maximum  likelihood  equation  is  [IS] 


P(g.T,r)  • L(g,T)  - r^(g  - I - e^)”^  + (100  - g - 


(78) 


- T 


- «.>■*] 


+ <T  ^ 

The  objective  function  for  the  generalised  three  •parameter  ganiia  die* 
trlbutlon  using  the  two  paranater  maximum  likelihood  equation  la  [IS] 


P(g^»r)  • l<(giar)  - t^(g  - 1 - «+)'^  + (100  - g - e^)"' 
+ (or  - 1-*^)  ^ + (100  • a • f^) 


(79) 
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rhe  objective  function  for  the  tvo-pareneter  game  distribution  using 


the  one  parameter  maximum  likelihood  equation  la  [IS] 

P(8,r)  - L(g)  - r[(g  - 1 - e^)"'  + (100  - g - e+)"^]  (80) 

The  algorithm  used  to  converge  the  objective  function  towards  the 
local  maximum  likelihood  is  as  follows  [IS]. 

1.  Maximise  the  objective  function,  P(‘T,r). 

2.  Check  for  convergence  to  the  optloam  l.e,  when 

|T(rj)  - T(r^_^)|  < e (81) 

where  < la  the  convergence  criterion  constant. 

3.  If  the  convergence  criterion  la  not  satisfied,  reduce  r^  by 
setting 

tj^.^  - drj,  0<  d<  1 (82) 

where  d la  a convergence  constant. 

4.  Increment  j and  repeat. 

The  maximisation  of  the  objective  function  has  been  done  by  many  non> 
linear  routines  [191,  However,  the  Hooka>Jeeves  pattern  search  method  [33] 
has  enjoyed  particularly  good  success  In  maximizing  MLB  objective  functions 
and  was  therefore  utilised  In  maximizing  the  objective  functions  for  the 
three - par ams ter  log  normal  distribution,  the  three •> par ams ter  Welbull  dis- 
tribution and  the  two,  three,  and  four- parameter  gaassa  distributions  [IS]. 

S.4  Goodness  of  Fit  Criteria 

Once  the  statistical  parameters  for  each  of  the  candidate  distribu- 
tions have  been  estimated,  the  distribution  which  the  data  follows  the 
closest  must  be  selected  from  the  candidate  distributions.  A statistical 
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method  which  Is  used  many  cliaes  to  find  out  how  well  data  fits  a certain 
distribution  Is  the  goodness  of  fit  test.  Several  goodness  of  fit  tests 
have  been  proposed  [121,  but  three  of  the  more  reliable  and  widely  used 
goodness  of  fit  tests  have  been  selected  as  criteria  for  the  selection  of 
the  "best"  distribution.  These  three  goodness  of  fit  tests  are  regres- 
sion, the  chi-square  test,  and  the  Kolmogorov- Smirnov  test. 

5. 4. a Regression 

Regression  In  its  simplest  form  involves  fitting  a polynomial  to  a 
set  of  given  data  plotted  on  certain  axes  [341.  In  the  case  of  fitting 
data  to  a statistical  distribution,  the  data  can  be  plotted  on  a plot 
whose  axes  correspond  to  certain  characteristics  of  that  particular 
statistical  distribution  (Section  5. 3. a).  It  Is  known  that  If  data  fol- 
lows that  particular  distribution,  then  the  data  will  follow  a straight 
line  fit  when  plotted  on  these  special  axes.  If  a linear  regression  Is 
performed  on  this  plotted  data,  It  can  be  determined  how  close  the  data 
does  fit  a straight  line.  This  then  provides  a measure  of  the  goodness 
of  fit  of  the  data  to  that  particular  distribution. 

If  a set  of  data  follows  the  two- parameter  normal  distribution,  a 
plot  of  the  data  with  the  X axis  as  a linear  scale  and  the  Y axis  as  a 
normal  probability  scale  will  follow  a straight  line  [18].  The  normal 
probability  scale  is  described  In  detail  In  Section  5. 3. a.  A typical 
plot  for  the  two- parameter  normal  distribution  Is  shown  In  Figure  11. 

If  a set  of  data  follows  the  two-parameter  log  noroial  distribution, 
a plot  of  the  data  with  the  X axis  as  a IoRj^q  ecale  and  the  Y axis  as  a 
normal  probability  scale  (Section  5. 3. a)  will  follow  a straight  line  [IS]. 
In  this  plot,  the  location  parameter  la  not  estimated  and  Is  assumed  to  be 
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IVt  5IHNDHRD  NORMAL  PROBABILITY  (PER  CENT) 


2-PflRflMETER 

NORMAL  DISTRIBUTION  PLOT 


Figure  11.  Typieel  2>P«r«a«e«r  Monuil  Dletrlbutloo  Plot 


2ero.  A typical  plot  for  the  two-parametar  log  normal  distribution  is 
shown  in  Figure  12. 

Both  the  plots  for  the  three>paramster  log  normal  distribution  and 
the  three-parameter  Weibull  distribution  have  been  discussed  in  Section 
5, 3. a.  A typical  plot  for  the  three-parameter  log  normal  distribution  is 
shown  in  Figure  13  and  a typical  plot  for  the  three -parameter  Weibull  dis- 
tribution is  shown  in  Figure  14.  The  three-parameter  gamna  distribution 
plot  requires  the  data  to  be  plotted  on  a plot  where  the  X axis  is  a 
linear  scale  and  the  Y axis  is  a gaosaa  probability  scale.  The  equation 
for  calculating  the  gamma  probability  scale,  H(s),  is  [27] 

F(x^)  r(g)  - t'^  e'*^  dt  (83) 

o 

where 

z - F(x^)  (84) 

where  F(x^)  is  the  cumulative  density  function  of  the  corrected  data  which 
is  given  by  equation  (58).  Equation  (83)  was  solved  iteratively  for  H(z) 
using  the  interval  halving  method  [271.  A typical  three -parameter  gasma 
distribution  plot  is  shown  in  Figure  15.  The  two-paramster  ganta  distri- 
bution plot  also  requires  the  X axis  to  be  a linear  scale  and  the  Y axis 
to  be  a gansna  probability  scale.  A typical  two-paraneter  gaoma  distribu- 
tion plot  is  shown  in  Figure  16.  In  each  of  the  above  plots,  the  data 
are  plotted  on  the  X axis  against  the  corresponding  median  ranks  on  the 
Y axis. 

Linear  regression  uses  linear  least  squares  which  uses  the  matrix 
approach  to  linear  regression  to  fit  a best  fit  straight  line  Co  Che  data 
[34].  Aa  a result  of  this  matrix  approach  to  linaar  regression,  a goodness 
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CUMULRTIVE  STANDARD  NORMAL  PROBABILITY  (PER  CENT! 


2-PflRflMETER 

LOG  NORMAL  DISTRIBUTION  PLOT 


N (CYCLES) 


Figure  12.  Typical  2*Paramecer  Log  Normal  Diatrlbution  Plot 
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CUMULATIVE  STRNOftRO  NORMAL  PROBABILITY  (PER  CENT) 


2-PflRflMETER 

LOG  NORMAL  DISTRIBUTION  PLOT 


N (CYCLES) 

Figure  12.  Typlcel  2*P«r«meeer  Log  Nornal  Dietrlbutlon  Plot 
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CUMULATIVE  STANDARD  NORMAL  PROBABILITY  (PER  CENT) 


3-PflRflMETER 

LOG  NORMAL  DISTRIBUTION  PLOT 


Figure  13.  Typlcel  3-FaraiMCer  Log  Normal  Dlatrlbutlon  Plot 
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WEIBULL  PROBABILITY  (PER  CENT) 


3-PflRflMETER 

WEIBULL  DISTRIBUTION  PLOT 


N ( CYCLES )-TBU  HAT 

Figure  lA.  Typical  3-FaraiMCer  Helbull  Distribution  Plot 
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GfiMMfl  PROBfiBILITY  (PER  CENT) 
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Figure  15.  Typical  3*Farain«Cer  GanBM  Distribution  Plot 
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GnMMft  PROBfiBILITY  (PER  CENT) 
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Figure  16.  Typical  2-ParaMCer  Gaona  Diacribution  Plot 
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of  fit  statistic,  called  the  coefficient  of  multiple  determination,  k , 

2 

can  be  calculated.  The  value  of  R is  always  between  zero  and  one.  The 

2 

closer  the  value  of  R is  to  one,  the  closer  the  fit  of  the  data  is  to  a 

2 

straight  line.  Therefore,  by  comparing  the  values  of  R for  each  of  the 

2 

distributions,  the  distribution  with  the  highest  value  of  R is  the  dis- 
tribution which  the  data  follows  the  closest. 

2 

This  value  of  R can  be  corrected  for  the  slope  of  the  least  squares 

line  in  an  attempt  to  achieve  a more  precise  measure  of  the  closeness  of 

2 

the  data  to  the  straight  line.  This  corrected  value  of  R is  called  the 

2 2 

closeness  and  is  given  the  symbol  C . The  derivation  of  C is  given  in 
Appendix  C. 


S.4.b  Chi-Square  Test 

The  chi-square  goodness  of  fit  test  is  a statistical  method  for  de- 
termining how  close  given  data  follow  a certain  distribution.  Basically, 


1)  divides  the  data  into  an  optimum  number  of  equiprobable  intervals, 

2)  counts  the  number  of  data  points  in  e.ich  interval  (called  the  ob- 
served frequencies), 

3)  calculates  the  number  of  data  points  that  should  be  in  each  in- 
terval based  on  the  estimated  distribution  parameters  (called  the 
expected  frequencies),  snd 

4)  compares  the  observed  frequencies  with  the  expected  frequencies 

[11,12]. 

2 

The  test  ststistlc,  x , * messurs  of  how  doss  the  observed  frequencies 

are  to  the  exi'ected  frequencies,  and  t'  us  how  close  the  data  follows  the 

2 

given  distribution,  x ia  given  by 
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^2-  I 

1-1  e, 

where  k is  the  number  of  equiprobable  intervals, 
are  the  observed  frequencies,  and 
e^  are  the  expected  frequencies. 

The  lower  the  value  of  the  chi-square  statistic,  the  closer  the  ob- 
served frequencies  match  the  expected  frequencies  and  thus  the  closer  the 
data  follows  the  given  distribution.  However,  the  chi-square  statistic 
can  not  be  compared  between  distributions  that  do  not  have  the  same  number 

of  distribution  parameters,  n.  because  the  degrees  of  freedom  for  the  chi- 

F 

square  statistic  for  distributions  not  having  the  same  number  of  distribu- 
tion parameters  Is  not  constant  [133.  Therefore,  the  tail  area  of  the 
chi-square  distribution  to  the  right  of  the  chi-square  statistic,  called 
A,  Is  computed  for  each  distribution  by  [13] 


f^2  exp(-u) 
A • ‘ 


(y-l) 

u du 
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where  v Is  the  number  of  degrees  of  frecdo::  and  u is  a variable  of  inte- 
gration. The  .slue  of  A Is  always  between  eero  and  one,  with  A equal  to 
"w  iciug  a perfect  fit.  The  lower  the  value  of  the  chl-'jquare  statistic, 
the  higher  the  value  of  the  tail  area,  all  other  thm.;  I'nr'  There- 
fore, the  distribution  lo  be  chosen  as  the  distribution  whi'.)  the  data 
follows  the  closest  is  th^  one  which  has  the  highest  value  of  A. 

The  chl-squara  statistic  may  be  ccrrparsd  with  a critical  value  which 
iollown  the  chi-square  distribution  at  an  acceptance  level  of  a with  v 
degrees  of  freedom,  x*  [12'',  .where 

a ^ sV 


Acceptance  of  the  proposed  distribution  as  the  distribution  which  the  data 
follows  should  occur  when  [12] 

2 2 

X i X (88) 
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The  tall  area,  A,  aay  be  compared  with  tha  acceptance  level  to  test  ac- 
ceptance of  the  proposed  distribution.  Acceptance  should  occur  when  [13] 

A a o (89) 

A 

The  end  points  for  the  classes  for  the  two  and  three-parameter  normal 
distributions  were  found  by  dividing  a standard  normal  curve  into  differ- 
ent numbers  of  equiprobable  intervals  [35].  The  end  points  for  the  equi- 
probable  intervals  for  the  three -paramo ter  Weibull  distribution  were 
given  by  [19] 

<1^  . f * i[.  l«{l  - (^}‘'']  (90) 

The  end  points  for  the  equiprobable  intervals  for  the  two,  three,  and 
four-parameter  gamma  distribution  were  given  by  [21] 

’i  ■ ^ ^ (i)}''*'] 

where  F.  ^ is  the  inverse  cumulative  density  function  for  the  generalited 
four-parameter  gamma  distribution. 

S.4,c  Kolmogorov-Smirnov  Test 

The  Kolmogorov- Smirnov  test  is  another  statistical  goodness  of  fit 
test  similar  to  the  chi-square  goodness  of  fit  test.  Basically,  it  cal- 
culates the  sample  cumulative  density  function  and  compares  it  with  the 
theoretical  cumulative  density  function  of  the  given  distribution  by  cal- 
culating the  maximum  deviation,  D,  between  the  two  cumulative  density 
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functions  [11].  The  test  statistic,  Z,  la  a measure  of  how  close  the  two 
cumulative  density  functions  are  and  thus  how  close  the  data  follows  the 
given  distribution  and  is  actually  equal  to  D. 

The  lower  the  value  of  the  Kolmogorov- Smirnov  statistic,  the  closer 
the  sample  cumulative  density  function  lies  to  the  theoretical  cumulative 
density  function,  and  thus  the  closer  the  data  follows  the  given  distribu- 
tion. Therefore,  the  distribution  to  be  chosen  as  the  distribution  which 
the  data  follows  the  closest  is  the  one  which  has  the  lowest  value  of  the 
Z statistic.  The  Kolmogorov- Smirnov  statistic  may  be  compared  with  a 
table  of  critical  values  to  determine  if  the  proposed  distribufibiT  Should 
be  accepted  as  the  distribution  which  the  data  follows  [36]. 
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SECTION  VI 


DETERMINATION  OF  THE  DISTRIBUTION 

Several  computer  programs  were  written  to  determine  the  distribution 
of  the  desired  fatigue  crack  propagation  variables  using  the  previously 
mentioned  statistical  concepts.  The  four  programs  written  to  determine 
statistical  distributions  of  fatigue  crack  propagation  variables  are: 

1}  Delta  N Distribution  Determination  Program  (Golden),  or 
DNDDPG, 

2)  Cycle  Count  Distribution  Determination  Program,  or  CCDDP, 

3)  Crack  Growth  Rate  Distribution  Determination  Program,  or 
CGRDDP,  and 

4)  Delta  N Distribution  Determination  Program  (HLE),  or  DNDDP. 

6. 1 Delta  N Distribution  Determination  Program  (Golden) 

This  program,  called  DNDDPG,  was  written  to  determine  the  distribu- 
tion of  the  4N/Aa  variable  computed  from  the  input  a vs.  N data  which  is 
supplied  by  program  DELTCP  (Section  7.1).  Basically,  it  fits  the  data  to 
four  distributions  and  computes  a goodness  of  fit  statistic  for  the  com- 
parison of  the  distributions.  The  four  distributions  fitted  are: 

1)  the  two-parameter  normal  distribution, 

2)  the  two-parameter  log  normal  distribution, 

3)  the  three -parameter  log  normal  distribution,  and 

4)  the  three-parameter  Wbibull  distribution. 
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It  uses  the  graphical  method,  including  the  Golden  Section  search  method, 

to  estimate  the  location  parameter  for  both  the  three- parama ter  log  normal 

distribution  and  the  three-parameter  Welbull  distribution.  The  goodness 

2 

of  fit  criterion  used  Is  C (Section  5. 4. a). 

This  program  produces  output  tfhlch  Includes  the  Input  a vs.  N data, 
the  computed  a vs.  N data,  some  of  the  test  conditions,  some  of  the  In- 
ternal program  parameters,  the  frequency  distribution  array,  the  dN/Aa 
data,  and  the  distribution  pai'ameters  and  a partial  analysis  of  variance 
table  for  each  distribution.  The  plots  generated  by  this  program  are  a 
relative  frequency  histogram,  a relative  cumulative  frequency  histogram, 
and  a distribution  plot  for  each  of  the  distributions.  Further  documenta- 
tion of  this  program  Is  shown  In  Appendix  D. 

6.2  Cycle  Count  Distribution  Deteiminatlon  Program 
This  program,  called  CCDDP,  was  written  to  determine  the  distribution 
of  the  N (cycle  count)  variable  from  a set  of  replicate  cycle  count  data 
at  one  crack  length  Level.  Identical  load  and  test  conditions  are  re- 
quired for  the  replicate  data.  This  program  fits  the  data  to  six  distri- 
butions. These  distributions  are: 

1)  the  two- parameter  normal  distribution, 

2)  the  two-parameter  log  normal  distribution, 

3)  the  three -parameter  log  normal  distribution, 

4)  the  three-parameter  Welbull  distribution, 

5)  the  three-parameter  gaona  distribution,  and 

6)  the  generalized  four- parameter  gamaa  distribution. 

It  uses  the  Maximum  Likelihood  Estimators  method  to  estimate  the  para- 
meters of  each  of  the  above  distributions  except  the  two-paraswter  normal 
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distribution  and  the  two-parameter  log  normal  distribution.  Three  good- 
ness of  fit  criteria  are  calculated  for  the  comparison  of  the  distribu- 
tions. They  are; 

1)  the  chi-square  tail  area, 

2)  Che  Kolmogorov- Smirnov  statlatic,  and 

2 

3)  R from  regression. 

This  program  produces  output  which  includes  the  input  replicate  cycle 
count  data,  the  test  conditions,  some  of  the  internal  program  parameters, 
the  frequency  distribution  array,  and  1)  the  estimated  distribution 
parameters,  2)  a partial  analysis  of  variance  table,  and  3)  the  goodness 
of  fit  criteria  for  each  distribution  except  the  generalised  four-para- 
meter  gamma  distribution,  for  which  only  the  estimated  distribution  para- 
meters and  the  goodness  of  fit  criteria  are  printed.  It  also  prints  a 
comparison  of  the  distributions  and  the  resulting  "best"  distribution 
based  on  the  goodness  of  fit  criteria.  The  plots  generated  by  this  pro- 
gram are  the  original  cycle  count  data  plot,  a relative  frequency  histo- 
gram, a relative  cumulative  frequency  histogram,  and  a distribution  plot 
for  each  of  Che  distributions  except  the  generalised  four-parameter  gamaa 
distribution.  Further  documentation  of  this  program  is  shown  in  Appendix 
E. 


6,3  Crack  Growth  Rate  Distribution  Determination  Program 
This  program,  called  CGRDDP,  was  written  to  determine  the  distribu- 
tion of  Che  crack  growth  rate  (da/dM)  variable  from  a set  of  replicate 
da/dN  data  at  one  crack  length  level.  This  da/dN  data  is  calculated  by 
Che  DADNCP  program  (Section  7.2).  Identical  load  and  Cast  conditions  are 
required  for  the  replicate  data. 
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This  program  is  nearly  identical  to  the  CCDDP  program  (Section  6.2), 
using  the  same  distributions,  the  same  parameter  estimation  method,  the 
same  goodness  of  fit  criteria,  and  having  nearly  the  same  output.  The 
main  difference  is  the  variable  of  interest  being  da/dN  instead  of  cycle 
count.  Thus  the  required  input  is  different  and  some  of  the  output  is 
different  in  this  respect.  Further  documentation  of  this  program  is 
shown  in  Appendix  F. 

6.4  Delta  N Distribution  Determination  Program  (MLg) 

This  program,  called  DNDDP,  was  vrittein  to  determine  the  distribution 
of  the  4N/Aa  variable  from  a sat  of  replicate  da/dN  data  et  one  crack 
length  level.  The  da/dN  data  used  is  the  same  as  that  used  by  the  CGRDDP 
program  (Section  6.3). 

This  program  is  based  on  the  CGRDDP  program.  One  main  difference  be* 
tween  them  is  that  the  input  da/dN  data  la  inverted  to  create  the  varl* 
able  AN/da.  The  second  main  difference  is  the  assumption  that  r for  the 
gaoma  distributions  is  equal  to  aero,  thus  reducing  the  3-parameter  gamma 
dittributlon  and  the  generalised  A-parameter  gamna  distribution  by  one 
parameter  (Section  5.2.]!).  Along  with  the  change  in  variable,  there  are 
appropriate  changes  in  the  output.  Further  documentation  of  this  program 
is  shown  in  Appendix  G. 
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SECTICM  VII 


GROWTH  RATE  AND  GROWTH  PREDICTION 


f 


Since  this  investigation  vas  not  Just  Interested  in  the  distribution 
of  fatigue  crack  propagation  variables  alone,  it  became  necessary  to 
write  several  other  programs  Co  aid  in  the  analysis  of  the  experitaenCal 
data.  These  supporting  programs  indued  1)  Delta  N Calculation  Program, 
or  DELTCP,  2)  da/dN  Calculation  Program,  or  DADNCP,  and  3)  a vs.  N Pre- 
diction Program,  or  AVNFRD.  Several  others  not  mentioned  here  were  used 
tc  in  the  analysi.e  and  manipulation  of  the  experimental  data. 

7.1  Delta  N Calculation  Program 

This  program,  called  DETJICP,  was  vrltten  to  calculate  Intermediate 
Aa  vs.  £hN  data  to  be  used  by  program  DMDDPG  (Section  6.1).  Basically, 
it  calculates  ^a  vs.  AN  data  from  a set  of  constant  amplitude  a vs.  N 
data  by  one  of  five  different  methods.  These  methods  are; 

1)  Che  secant  method, 

2)  reject  certain  selectable  data  points  and  use  the  secant  method, 
thereby  increasing  As, 

3)  the  quadratic  7-poinc  incremental  polynomial  method, 

4)  reject  certain  selectable  data  points,  recreate  new  a vs.  N data, 
and  then  use  the  quadratic  7-point  incremantsl  polynomial  method, 
and 

5)  use  Che  quadratic  7-point  incremental  polynomial  method,  recreate 
new  a vs.  N data,  reject  certain  selectable  data  points,  and  than 
use  the  secant  method. 


4 
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Further  documentetion  of  this  program  Is  shown  in  Appendix  H. 

7 .2  da/dN  Calculation  Program 

This  program,  called  DADNCP,  was  written  to  calculate  the  crack 
growth  rate,  da/dN,  by  the  six  different  methods  presented  in  Section 
4.  These  methods  are; 

1)  the  secant  method, 

2)  the  modified  secant  method, 

3)  the  linear  7*poinC  incremental  polynomial  method, 

4)  the  quadratic  7-point  incremental  polynomial  method, 

5)  the  linear  log-log  7-point  incremental  polynomial  method,  and 

6)  the  quadratic  log- log  7-poinc  Incremental  polynomial  method. 

For  each  of  these  methods,  the  calculated  da/dN  data  is  integrated  back 
into  estimated  a vs.  N data,  which  is  compared  with  the  original  a vs.  N 
data,  resulting  in  an  average  incremental  error.  By  comparing  these 
errors,  the  da/dN  calculation  method  which  results  in  the  lowest  error 
can  be  selected. 

The  required  input  for  this  program  is  a set  of  constant  Aa  a vs.  N 
data.  This  program  produces  output  which  includes  the  input  a vs.  N data, 
the  test  conditions,  da/dN  vs.  AK  and  actual  cycle  count  data  vs.  esti- 
mated cycle  count  data  for  each  da/dN  calculation  method,  and  a sunoary 
of  the  errors  from  each  method  with  the  resulting  "best"  da/dN  calculation 
method.  Further  documentation  of  this  program  is  shown  in  Appendix  I. 

7.3  a vs.  N Prediction  Program 

This  program,  called  AVNFRD,  predicts  a vs.  N data  from  the  distri- 
bution of  da/dN  (or  dN/da)  as  a function  of  crack  length  and  compares  it 


60 


with  ch«  original  a va.  N data.  Tha  raqulrad  input  ia  tha  knowladga  «( 
tha  dlatrlbutlon  of  da/dN  (or  dft/da)  aa  a funation  of  araak  langck  aa 
daearmlnad  by  tha  COkODf  (or  DNDDf)  progran.  Thii  profraai  aateaCi  • 
growth  rata  at  aach  craak  langth  ualng  a randaa  ganaraior  and  tha 

dlatrlbutlon  paramtari.  Thla  growth  rata  la  ttian  uaad  to  aalmiUta  gN 
aa  a function  of  craek  langth,  which  ic  uaad  to  pradlce  rapllaata  eata 
of  a va.  N data.  Thaaa  pradlatad  aata  of  a va.  II  data  a«*a  than  aoiaparad 
with  tha  original  a va.  M data  auta. 

Thla  prograai  produeaa  output  which  Includao  tha  tact  conditional  tha 
pradletad  a va.  M data,  and  a plot  of  all  of  tha  pradlatad  a va.  N datti 
further  docuaiantatlen  of  thla  prograa  la  ahown  In  Appundla  J, 


II 


with  the  original  a vs.  N data.  The  required  input  is  the  knov ledge  of 
the  distribution  of  da/dN  (or  dN/da)  as  a function  of  crack  length  as 
determined  by  the  CGRDDP  (or  DNDDF)  program.  This  program  selects  a 
growth  rate  at  each  crack  length  using  a random  number  generator  and  the 
distribution  parameters.  This  growth  rate  is  then  used  to  calculate  AN 
as  a function  of  crack  length,  which  is  used  to  predict  replicate  sets 
of  a vs.  N data.  These  predicted  sets  of  a vs.  N data  are  then  compared 
with  the  original  a vs.  N data  sets. 

This  program  produces  output  which  includes  the  test  conditions,  the 
predicted  a vs.  N data,  and  a plot  of  all  of  the  predicted  a vs.  N data. 
Further  documentation  of  this  program  is  shown  in  Appendix  J. 


SECTION  VIII 


STATISTICAL  ANALYSIS  OF  PREVIOUSLY  GENERATED  DATA 

A considerable  amount  of  crack  propagation  data  In  the  form  of  a ve. 
N data  have  recently  been  generated  at  Purdue  University  for  center  crack 
specimens  of  2024-T3  aluminum  alloy  [37].  From  this  set  of  data,  there 
were  30  different  over load /under load  tests  which  were  conducted  under 
constant  stress  Intensity  conditions  and  at  constant  Aa.  From  each  of 
these  tests,  approximately  19  to  155  data  points,  for  a total  of  2076 
data  points,  were  collected  after  the  crack  had  grown  through  the  region 
influenced  by  the  over  load /under load  sequence.  The  data  typically  chosen 
for  analysis  Is  shown  In  Figure  17.  This  large  amount  of  data  was  col- 
lected following  the  overload  affected  region  to  establish  a final  steady 
state  growth  rate  as  well  as  establishing  the  steady  state  growth  rate 
for  the  next  test  [37,38,39].  From  this  set  of  test  results,  there  are 
2 to  7 sets  of  data  at  each  of  five  different  loading  conditions. 

The  value  cf  these  data  for  statistical  evaluation  centers  on  the 
accuracy  with  which  the  original  a vs.  N data  were  collected.  In  these 
tests,  the  crack  length  was  monitored  and  measured  with  a lOOX  micro- 
scope mounted  on  a digital  measureoenc  traverse.  The  traverse  has  a 
resolution  of  0.001  am  (0.00004  In.)  with  a direct  digital  read-out.  A 
printer  activated  by  a push  button  was  connected  to  the  cycle  counter 
and  the  digital  traverse.  In  collecting  the  data,  the  microscope  was 
advanced  an  Increment  of  0.01  mn,  0.02  an,  or  0.05  mm  (depending  on  the 
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:ngth  in  h 


growth  rate).  When  thf>  crack  had  grown  this  Incrcnent  as  observed  with 
the  cross  hair  In  the  Er.lcroscope , the  printer  was  activated  with  the  push 
b'ltton  and  the  crack  length  and  number  of  cycles  were  printed.  The  re- 
sulting data  are  very  dense  and  appear  to  be  fairly  accurate.  This  large 
amount  of  data  was  used  to  make  a preliminary  statistical  analysis  to  aid 
in  the  direction  and  scope  of  this  lnv>  >.itlon  [40]. 

8 . 1 Distribution  of 

The  first  step  of  the  analysis  was  to  determine  the  distribution  of 
the  variable  ^N/da  which  was  calculated  by  the  secant  method.  This  was 
done  by  writing  a pair  of  programs  using  many  of  the  statistical  concepts 
presented  In  Section  5.  These  programs,  called  Delta  N Calculation  Pro- 
gram, or  DELTCF  (Section  7.1),  and  Delta  N.  Distribution  Determination 
Program  (Golden),  or  DNDDPG  (Section  6.1),  were  run  on  each  of  the  data 

sets.  The  distributions  were  ranked  from  1 to  4 (1  being  the  best)  based 

2 

on  the  goodness  of  fit  criterion,  C (Section  5. 4. a).  The  rankings  were 
averaged  over  all  of  the  tests  and  the  results  are  shown  in  Table  1.  The 
best  distribution  was  the  three-parameter  log  normal  distribution  fol- 
lowed closely  by  the  two-parameter  log  normal  distribution.  A plot  cf 
the  fit  of  the  dN/de  data  to  the  three -parameter  log  normal  distribution 
is  shown  In  Figure  18. 

Based  on  these  results  and  the  use  of  the  DBLTCP  and  DHTOPC  programs, 
the  follcMing  conclusions  were  sMde. 

1)  Ilie  2-paramster  Welbull  distribution  was  tried  and  re- 
l<3cted  fi<M  all  further  analysis  because  of  its  poor  per- 
forman.;  ; providing  a fit  for  the  dN/Aa  data  due  to  it's 
..ack  of  j location  parameter. 
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2}  Include  the  other  four  dlst'-ibutlons  in  the  analysis  of 
other  fatigue  crack  propagation  variables. 

3)  Using  the  constant  amplitude  portion  of  overload/underload 
data  does  not  lead  to  a satisfactory  statistical  analysis. 
Therefore,  a statistically  designed  test  program  was  needed. 

A)  The  graphical  method  of  parameter  estimation  tended  to  be 
unstable  and  unreliable  for  the  data  used.  Therefore,  the 
Maximum  Likelihood  EstlmaLc-rs  method  of  parameter  estlma* 

tlon  was  tried  and  used. 

2 

5)  The  use  of  C as  a goodness  of  fit  criterion  was  poor  be- 
cause  It  failed  to  distinguish  between  the  distributions 
very  well.  Therefore,  the  chl*squere  ar  iiolnogorov*Smirpov 
goodness  of  fit  tests  were  tried  and  used. 

8.2  Effect  of  Quadratic  7-Polnt  Incremental  Polynomial  Mstltod 

/ ‘ 

The  second  step  of  the  analysis  was  to  examine  the  effet  using 
Che  quadratic  7>polnt  Incremental  polynomial  method  vs.  using  t recant 
method  In  calculating  the  variable  dN/Aa.  This  was  done  by  runiu..*  the 
DELTCP  program  and  changing  the  calculation  method  for  each  of  the 
data  aacs.  Once  the  AN/da  data  was  calculatad  for  each  data  eet,  i.c  was 
run  on  the  DNDDPC  program  Co  datarmlna  tha  affact  of  the  AN  calcuLi  on 
method  on  the  distribution  pareassters.  The  most  noticeable  effect  wee 
the  decrease  In  Che  variance  using  the  quadratic  7'polnt  incremental 
polynomial  machod  as  shown  in  Tabla  II.  From  this,  it  is  avidant  that 
tha  quadratic  7«point  Incramsntal  polynomial  mathod  Introduces  quite  a 
savothlng  affect  In  reducing  Che  amount  of  data  scactar  and  thus  the  data 
varlanca . 
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Xabla  II 


Smoothing  Effect  o£  the  Incremental  Polynoalel  Method 


I 

DISTRIBUTION 


nvomsz 


VfiR. 

vm.  (secrktT 


STD. 


DEV.  OF 


vm.  (i.r.) 

vm.  (SECflNT) 


2-PflRflMETER 

NORMRL 

0.H19 

U.0820 

2-PflRflMETER 
LOG  NORMAL 

0.H44 

0.1273 

S-PflRaHETER 
LOG  NORMAL 

0.647 

0.2817 

3-PARAMETER 
WE I BULL 

0.871 

0.4185 

8.3  Life  Prediction  Using  Estimated  Dtetrlbution  Parameters 


The  next  step  In  the  analysis  was  to  see  if  the  estimated  distri- 
bution parameters  could  be  used  for  life  prediction.  Using  the  mean  c£ 
the  ANMa  data  (for  the  two-parameter  normal  distribution)  and  the  over- 
all change  In  crack  length  (a^-a^) , the  final  cycle  count,  N^,  was  pre- 
dicted and  compared  with  the  observed  value  of  for  each  set  of  data 
and  then  averaged  over  all  the  data  sets.  The  results  are  shown  in  Table 
III.  From  the  lalatively  low  amount  of  error,  It  is  evident  that  statis- 
tical methods  using  estimated  distribution  parameters  could  prove  Invalu- 
able for  life  prediction. 


Table  III 

Life  Prediction  Based  on  th«  Mean 


N, 


flVERRGE  ^ 


PERCENT  ERROR 


8.4  Effect  of  Aa 

The  final  step  In  tho  analysis  of  toe  previously  generated  data  was 
to  determine  the  affect  of  the  else  of  This  was  done  by  using  Che 

DELTCF  program  to  generate  AN  data  with  different  values  of  Aa.  By 
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rejecting  certain  successive  data  points  (i.e.  every  1 out  of  2,  every  2 
out  of  3,  etc.)»  data  with  Increasing  values  of  Aa  were  generated.  The 
DNDDPG  program  was  then  run  on  each  different  ^a  set  of  data  for  each  of 
the  data  sets.  Also,  several  tests  at  the  sasie  load  conditions  were  com- 
bined to  give  a large  amount  of  data  and  then  A a was  Increased  as  de- 
scribed above.  The  results  are  shown  in  Figure  19  and  Table  IV.  From 
these  results,  it  is  obvious  that  the  larger  As  is,  the  smaller  the  re- 
sulting variance  of  the  data  will  be. 
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Table  IV 

Effect  of  Increaalng  La 


VW.  (^.10  ifl)  VW.  (fiife.18  m) 

DISTRIBUTION 


2-PflRflMETER 

NORMAL 

0.617 

0.513 

2-PflRflHETER 
LOG  NORMAL 

0.660 

0.S39 

3-PflRflMETER 
LOG  NORMAL 

0.851 

0.412 

3-PflRflMETER 
WE I BULL 


0.883 


0.761 


SBCTIOK  IX 


BmRi>eNrAL  investigation 

In  an  effort  to  answer  the  investigation  objectives,  it  becaoie 
necessary  to  conduct  an  experiaental  investigation  to  provide  adequate 
data  for  subsequent  analysis.  Through  the  use  of  previously  collected 
data  (Section  8),  it  becaaie  increasingly  clear  that  any  experiiaental  in- 
vestigation that  would  be  expected  to  provide  neanlngful  results  would 
have  to  be  statistically  designed.  Through  the  use  of  sone  preliminary 
theoretical  and  experimental  testing,  a test  program  was  designed. 

9. 1 Experimental  Test  Proaram 

Given  the  objectives  of  the  investigation  (Section  3),  it  was  evi- 
dent that  replicate  tests  under  Identical  load  and  environmental  condi- 
tions had  to  be  conducted.  It  was  also  obvious  that  constant  asiplltude 
loading  should  be  used  rather  than  constant  AK  (load  shed)  loading  since 
it  would  be  much  easier  to  control  and  replicate  and  also  give  a range  of 

AK  levels.  To  be  able  to  find  the  distributions  of  N and  da/dN,  the  data 

from  each  test  had  to  be  taken  at  consistent  discrete  a levels. 

To  detersiine  the  actual  load  levels  to  be  used,  several  preliminary 
tests  using  the  sane  lot  of  the  sane  naterial  ware  conducted.  To  obtain 
the  desired  growth  rates  ^ x lO'^  in. /cycle  and 

da/dNf^ajf  S!  5 X 10~^  in. /cycle)  and  keep  the  testing  tine  within  reason, 
it  was  found  that  AP  should  bo  4200  lbs.  It  was  also  determined  to  use 

an  R ratio  of  0.2  to  stay  well  out  of  the  compression  region. 
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A preliminary  theoretical  investigation  was  conducted  to  determine 
where  the  data  was  to  be  taken.  It  was  found  that  to  get  the  desired 
range  of  growth  rates,  the  data  would  have  to  be  taken  over  at  least 
40.0  nm.  It  was  determined  chat  steady  state  conditions  would  not  exist 
until  9.0  ma  due  to  the  crack  initiation  load  shedding  process.  In  an 
effort  to  reduce  data  error  as  much  as  possible  and  still  obtain  a rea- 
sonable amount  of  data, the  initial  Aa  was  chosen  to  be  0.20  nm  based  on 
the  statistical  analysis  of  previous  data  (Section  8.4).  Since  the 
growth  rate  would  be  Coo  fast  to  operate  Che  optical  system  and  the 
printer  at  the  end  of  the  test  for  the  load  levels  chosen,  Aa  would  be 
increased  to  0.40  nm  and  finally  to  0.80  nm.  The  number  of  data  points 
taken  at  Aa  0.40  nm  and  As  • 0.80  nm  mre  arranged  so  that  when  succes- 
sive data  points  were  rejected  (to  find  the  effect  of  increasing  As), 
there  would  be  no  large  gaps  in  the  data.  A schematic  representation  of 
the  test  program  is  shown  in  Figure  20. 

In  order  to  obtain  enough  data  to  conduct  a meaningful  statistical 
analysis,  it  was  determined  that  there  should  be  at  least  50  replicate 
tests  [13].  However,  since  more  specimens  were  available,  a total  of  68 
tests  were  conducted,  thereby  Increasing  the  confidence  of  the  statistical 
analysis  results.  The  test  conditions  are  listed  below. 

Sq  ■ 9.00  am. 
a^  “ 49.80  mm. 


R • 0.20 

P , - 1050  lbs. 

min 

P - 5250  lbs. 
max 

AP  - 4200  lbs. 
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9.2  Test  Spectnen 

The  test  specimens  used  in  this  investigation  were  O.lOO  inch  thick 
center  crack  panels  of  2024-T3  aluminum  alloy.  The  specimen  geometry  is 
shown  in  Figure  21. 

Test  specimens  were  obtained  with  a mill  finish  and  polished  to  a 
mirror  finish  In  the  vicinity  of  the  crack  path  to  facilitate  optical 
observation  of  the  crack  tip  during  crack  growth  measurement.  The  lot  of 
specimens  was  numbered  in  order  as  they  were  taken  out  of  the  shipping 
crate  so  that  true  randomisation  of  the  samples  could  be  accomplished. 

The  fixture  plate  holes  were  drilled  and  reamed  to  the  desired  dimen- 
sions.  The  stress  raiser  shown  in  detail  In  Figure  2i  was  machined  with 

I 

an  electro-discharge  machine. 

Before  loading  each  specimen,  the  centerline  of  the  specimen  was 

t 

scribed  at  the  stress  raiser  and  a silica  gel  desiccant  was  applied  at 
»he  stress  raiser.  The  entire  expected  crack  path  was  then  sealed  with 
clear  polyethylene  to  insure  desiccated  air  at  the  crack  tip.  Loading 
was  then  applied  parallel  to  the  direction  of  rolling  of  the  material. 


I 
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9.3  Test  Equipment 

The  test  machine  was  a 20  Kip  electro-hydraulic  closed- loop  system 
operated  In  load  control.  A function  generator  waa  used  to  generate  a 
sinusoidal  voltage  signal  which,  when  superls^osed  on  a d.c.  set  point 
voltage,  constituted  the  desired  input  to  the  system.  During  testing,  sn 
oscilloscope  was  used  co  mucltor  the  feedback  signal  (load)  and  tba  output 
of  Che  aaplicwte  isnA.suremer.c  system  of  the  testing  machine  to  insure  cor- 
rect load  lavals  Mjd  sinusoidal  loading.  A digital  cycle  counter  was  used 
to  count  the  nuoiber  of  applied  load  eye  Us.  Crack  growth  was  monitored 
with  a zoom  stereo  microscope  operated  at  a magnification  of  ISOx 
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rigidly  mounted  on  a horieontal  and  vertical  digital  traversing  system. 

A crosshair  mounted  in  the  microscope  was  used  as  a reference  line  during 
data  acquisition.  A digital  resolver  system  on  the  horiEontal  traverse 
produced  a digital  output  vith  a resolution  of  0.001  nm  (.00004  In.}. 

The  direction  of  travel  of  the  optical  system  prior  to  data  acquisition 
was  never  changed  during  a test  to  eliminate  any  hysteresis  effects  in 
the  traverse  system.  Both  the  digital  traverse  and  cycle  counter  outputs 
(crack  length  and  number  of  cycles)  were  connected  to  a mechanical 
printer.  The  printer  printed  both  the  crack  length  and  the  cumulative 
cycle  count  by  the  operation  of  a push  button.  A strobe  light  synchro- 
nised with  the  feedback  signal  was  triggered  at  the  point  in  the  load 
cycle  when  the  crack  was  most  fully  open  to  illuminate  the  crack  tip. 

More  detailed  discussions  of  the  test  equipoent  can  be  found  in  refer- 
ences [37,38,39]. 


9.4  Test  Procedure 

Since  the  scope  of  this  investigation  strictly  involved  the  deter- 
mination of  the  effect  of  material  properties  on  fatigue  crack  propaga- 
tion, care  was  taken  to  control  as  many  other  variables  as  possible.  All 
tests  were  subject  to  nearly  identical  environmental  conditions  of  room 
temperature  (24°C)  and  desiccated  air.  Loads  wre  controlled  to  within 
0.27.  of  the  desired  load  using  the  test  machine's  aa^litude  maasuremant 
system.  To  prevent  any  effects  from  the  order  in  which  the  specimens 
were  run,  the  specimens  were  randomised  using  a computer  program  which 
utilized  a random  number  generator.  The  tests  were  run  in  the  random 
order  determined  by  this  program.  The  order  of  teste  is  shown  in  Appen- 
dix K. 
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Crack  Initiation  starting  at  the  stress  raiser  was  performed 
starting  at  ■ 15000  lbs.  and  shedding  the  load  107.  no  sooner  than 
every  0.5  no  (12.5  times  the  change  in  plastic  zone  radius  due  to  the 
load  shed)  to  the  desired  test  load  level.  Fatigue  cycling  was  done 
initially  at  10  hz  up  to  5.4  ora  (due  to  reduced  frequency  response  of 
the  testing  machine  at  high  loads)  and  then  at  20  hz.  To  make  certain 
that  no  load  effects  were  present  in  the  data,  the  test  load  level  was 
reached  1.0  ora  before  data  acquisition  (58  times  the  change  in  plastic 
zone  radius  due  to  the  last  load  shed).  The  load  level  was  held  constant 
throughout  the  test  (thus  increasing  AK  with  increasing  crack  length). 

All  tests  were  started  at  the  same  init'a!  crack  length  (2a  ■ 18.00  ran). 
The  location  of  the  centerline  of  the  specimen  was  noted  as  a reference 
to  insure  consistent  crack  length  measurements  throughout  the  test. 

Cycling  was  continuous  throughout  the  test  to  eliminate  any  time  or 
underload  effects  on  subsequent  fatigue  crack  growth. 

The  crack  length  and  number  of  cycles  were  monitored  continuously 
for  each  test  and  discrete  data  points  were  taken  as  determined  by  the 
test  program.  Data  were  actually  taken  by  advancing  the  optical  system 
by  Che  specified  increment  and  pressing  Che  printer  push  button  when  the 
crack  tip  had  grown  to  Che  incremented  position  as  determined  by  the 
crosshair  in  the  stereo  microscope.  The  amount  of  error  in  Che  data  ac* 
quisicion  process  is  given  in  Section  9.5. 

9.5  Metosuiewent  Accuracy 

In  an  attempt  to  Isolate  Che  data  variance  due  to  the  material  prop- 
erties, a measure  of  the  experimental  error  was  needed.  This  axperimental 
error  results  from  the  random  error  in  measuring  the  cycle  count  and  the 
crack  length. 
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By  using  the  test  machine's  anplltude  neasuremant  system  which  com- 
pares a known  input  signal  with  the  feedback  signal  (applied  load),  the 
loads  can  be  controlled  to  within  0.2%. 

Brror  in  the  crack  length  measurement  is  due  to  two  sources.  If  the 
spatial  relationship  between  the  microscope  crosshair  and  the  scribed 
reference  line  on  the  specimen  is  not  constant,  then  an  undetermined 
amount  of  measurement  error  is  present.  This  usually  occurs  when  the 
microscope  is  accidentally  moved  with  respect  to  the  specimen  and  can  be 
avoided  by  a careful  experimental  procedure. 

The  second  source  of  crack  length  measurement  error  is  the  alignment 
of  the  crack  tip  with  the  microscope  crosshair.  This  alignment  process 
consists  of  1)  defining  the  crack  tip  location,  2)  defining  the  crosshair 
location,  and  3^  comparison  of  the  two  locations  to  see  if  they  are  iden- 
tical. If  they  are,  then  the  printer  button  is  pushed  and  a data  point 
is  taken. 

To  determine  how  well  the  observer's  eye  performs  this  alignment  pro- 
cess, the  following  test  was  devised.  A crack  was  initiated  and  the  cy- 
cling was  stopped  when  the  observer  determined  that  the  crack  had  reached 
9.00  mm.  He  then  took  10  repeat  measurements  of  the  crack  length,  being 
careful  to  always  approach  the  crack  tip  from  the  same  direction  to  pre- 
vent any  hysteresis  effects.  This  series  of  10  repeat  measurements  was 
repeated  at  9 other  predetermined  crack  lengths.  The  mean  and  standard 
deviation  of  each  set  of  10  repeat  measurements  was  computed  and  the 
error  of  the  original  data  point  was  then  calculated  in  terms  of  the 
standard  deviation.  The  results  of  the  10  seta  of  repeat  meeaurements 
are  as  follows. 
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• 0.001414  mn. 

Sg  • 0.001390  nra. 

where 

Xg  Is  the  mean  of  the  errors, 

is  the  standard  deviation  of  the  errors. 

Therefore,  the  average  experimental  error  for  each  data  point  is  0.001414 
mm.  The  average  experimental  error  as  a function  of  the  crack  length 
measurement  interval,  4a,  is  shown  in  Table  V.  It  should  be  noted  here 
that  the  larger  Aa  is,  the  smaller  the  average  experimental  error  is. 


i 
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Table  V 


Average  Experimental  Error 


Afl  INCREMENT  (MM)  flVERflSE  ERROR  (PERCENT) 


0.20 

0.71 

0.40 

0.35 

0.80 

1 

0,17 

SECTION  X 


DATA  ANALYSIS  AND  RESULTS 

As  a result  of  the  experimental  Investigation  conducted  as  described 
In  Section  9,  68  replicate  a vs.  N data  sets  were  obtained.  These  data 
are  shown  in  Figure  22.  Using  these  data,  an  analysis  was  performed  to 
meet  the  objectives  of  the  investigation  (Section  3). 

10.1  Distribution  of  N 

The  first  objective  to  be  net  was  to  determine  the  distribution  of  N 
as  a function  of  crack  length.  The  replicate  N data  used  was  readily  ob- 
tained from  the  original  replicate  a vs.  N data.  Typical  replicate  cycle 
count  data  are  shown  in  Figure  23.  The  distribution  of  the  replicate 
cycle  count  data  was  determined  at  each  crack  length  level  through  the 
use  of  Che  CCDDP  program  (Section  6.2).  At  each  crack  length  level,  this 
program  calculated  Che  distribution  parameters  and  goodness  of  fit  criteria 
for  Che  six  distributions  and  then  compared  the  goodness  of  fit  criteria 
between  five  of  the  distributions  in  order  Co  establish  Che  distribution 
rankings.  The  generalized  A-parameter  gamns  distribution  was  not  consid- 
ered for  the  distribution  rankings  because  it  wss  expected  to  have  an  ex- 
cellent fit  to  the  cycle  count  data  due  to  it's  power  parameter  (Section 
S.2.e).  The  distribution  parameters,  goodness  of  fit  criteria,  and  the 
distribution  rankings  ware  then  combined  over  all  of  Che  crack  length 
levels . 
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Figure  22.  Orlgloel  Replicate  a ve.  M Data 


84 


fl  (CRRCK  LENGTH  IN  MM) 


37.00 


36.80 


fl  VS.  N 

REPLICATE  Cfl  TESTS. 

DELTA  P = 4.20  KIP 
68  DATA  POINTS 
AO  r 9.000  MM 
R = .20 


36.60H 


36.H0H 


36.20H 


XOC  XX  X X 


X 


36.00 H 


3S.80-I 1 1 1 1 1 

.2000  .2200  .2H00  .2600  ^ .2800  .3000 

N (CYCLES)  (XIO  ®) 


Figure  23.  Typical  Replicate  Cycle  Count  Data 
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The  distribution  parameters  of  the  cycle  count  data  as  a function  of 
crack  length  were  plotted  for  each  of  the  six  distributions  and  are  shown 
in  Figures  2^  through  29.  The  distribution  parameters  are  normalised  so 
that  their  minimum  and  maximum  values  are  equal  to  zero  and  one,  re* 
spectlvely.  As  a result  of  this  normalization,  these  figures  do  not  show 
the  actual  values  of  the  distribution  parameters  but  are  intended  to  re- 
flect trends  present  in  these  parameters. 

The  goodness  of  fit  criteria  for  each  distribution  wern  averaged  over 

all  of  the  crack  length  levels.  These  results  are  shown  in  Table  VI.  For 

these  goodness  of  fit  criteria,  the  best  fit  of  the  data  to  a distribution 

occurs  when  the  chi*square  tail  area  is  a maximum,  the  Kolmogorov-Smirnov 

2 

statistic  is  a minimum,  and  the  closeness,  R , is  a maximum.  Using  these 
relationships,  an  understanding  of  which  distributions  provide  the  best 
fit  for  the  cycle  count  data  can  be  obtained. 

Tne  distribution  rankings  at  each  crack  length  level  were  combined 
over  all  of  the  crack  length  levels.  By  convention,  the  lower  the  value 
of  the  distribution  ranking,  the  better  the  fit  of  the  data  to  the  given 
distribution.  The  mean  rank  and  it’s  standard  deviation  for  each  of  the 
distributions  and  the  number  of  times  each  distribution  was  selected  as 
the  best  distribution  were  calculated  during  this  combining  process. 

These  results  are  shown  in  Table  VII. 

The  3-parameter  log  normal  distribution  provided  the  best  fit  for  the 
cycle  count  data  by  a wide  margin,  as  evidenced  by  Che  low  distribution 
ranking  value,  the  low  Kolmogorov-Smirnov  test  statistic  value,  and  the 
very  large  number  of  times  it  was  selected  as  the  best  distribution.  The 
3-paraffleter  gamma  distribution  provided  the  next  best  fit,  while  the  2- 
pararoeter  log  normal  distribution  and  the  3-paraineter  Weibull  distribution 
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Figure  24.  Z-ParecpeCer  Norokal  Distribution  Parameters 
of  Cycle  Count  Data  as  a Function  of  Crack 
Length 
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o£  Cycle  Count  Data  as  a Function  of  Crack 
Length 
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NORMfiLHED  DISTRIBUTION  PflRflHETERS 


Table  VI 

Average  Goodness  of  Fit  Criteria  for  the 
Distribution  of  Cycle  Count  Data 


CHI-SQURRE  KOLMOGOROV-  CLOSENESS 

DISTRIBUTION  tail  rrer  smirnov  test  ir  squared) 


0.93310 


2-PRRRMETER 

NORMRL 

0.8365 

0.0995 

2-PflRflMETER 
LOG  NORMAL 

0.8842 

0.0857 

3-PflRflMETER 
LOG  NORMAL 

0.8594 

0.0699 

3-PFIRflMETER 

WE I BULL 

0.8340 

0.0882 

3-PflRflMETER 

GRMMR 

0.8602 

0.0722 

GENERALIZED 

4-PRRRMETER 

0.8075 

0.0722 

GRMMR 

0.95799 


0.98223 


0.93658 


0.97160 


Table  VII 

Distribution  Rankings  for  the  Distribution  of  Cycle  Count  Data 


DISTRIBUTION  mean 


NUMBER  OF 
STANDARD  TIMES  BEST 
DEVIATION  DISTRIBUTION 


2-Pf)RflMETER 

NORMAL 

4.982 

0.1348 

a-PRRRMETER 
LOG  NORMAL 

3.147 

0.6780 

3-PflRflMETER 
LOG  NORMAL 

1.221 

0.5882 

3-PRRflMETER 
WE I BULL 

3.650 

0.6338 

3-PHRRMETER 

GAMMA 

2.000 

0.4969 

tied  for  the  third  best  fit  for  the  data.  The  2-paraiiieter  normal  distri- 
bution finished  last  in  the  distribution  ranking  as  it  provided  a very 
poor  fit  for  the  data. 

10.2  Crack  Growth  Rate  Calculation  Methods 

The  second  objective  to  be  nt  was  to  determine  which  crack  growth 
rate  calculation  isethod  introduced  the  least  amount  of  error  into  Che 
da/dN  data.  This  was  to  be  done  by  integrating  the  da/dN  data  calculated 
by  each  crack  growth  rate  calculation  method  back  into  a vs.  N data  and 
Chen  calculating  the  error  between  the  new  a vs.  N data  and  the  original 
a va.  N data. 

The  DADNCP  program  (Section  7.2)  was  run  on  each  of  Che  68  original 
a vs.  M data  sets.  This  program  calculates  the  da/dN  vs.  dK  data,  inte** 
grates  the  da/dN  data  back  into  a vs.  N data  using  Simpson's  one-third 
rule  and  the  Crapeeoidal  rule,  and  Chen  calculates  a step  by  step  average 
incremental  error,  as  outlined  by  Frank  and  Fisher  [2],  for  each  of  the 
six  da/dN  calculation  methods.  The  da/dN  calculation  method  which  re- 
sults in  the  lowest  average  incremental  error  is  then  selected  as  the 
best  da/dN  calculation  method  for  that  data  set.  The  log|^Q  da/dN  vs. 
log^O  data  are  plotted  for  each  of  the  da/dN  calculation  methods  and 
typical  plots  of  these  data  are  shown  in  Figures  30  through  33. 

The  average  incresMntal  error  from  each  da/dN  calculation  method  was 
averaged  over  all  of  cha  data  sets  and  the  number  of  times  each  da/dN 
calculation  mschod  was  selected  es  the  best  msthod  wes  computed.  These 
results  ere  shown  in  Teble  Vlll.  The  modified  secent  msthod  hed  the  lowest 
average  incremental  error,  followed  closely  by  the  secant  swthod.  The 
modified  secent  method  end  the  secent  msthod  were  both  selected  es  best 
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table  VllI 


de/dN  Calculation  Method  Result s 


Dfl/DN 

CRLCULRTION 

METHOD 

OVERALL  AVERAGE 
INCREMENTAL  ERROR 
(PER  CENT) 

NUMBER  OF 
TIMES  BEST 
METHOD 

SECANT  METHOD 

2.70 

17 

MODIFIED 

SECANT  METHOD 

2.58 

51 

LINEAR 

7-PaiNT  INCREMENTAL 
POLYNOMIAL  METHOD 

6.96 

0 

QUADRATIC 

7-POINT  INCREMENTAL 
POLYNOMIAL  METHOD 

6.83 

0 

LINEAR  LOG-LOG 
7-POINT  INCREMENTAL 
POLYNOMIAL  METHOD 

• 

CD 

0 

QUADRATIC  LOG-LOG 
7-PaiNT  INCREMENTAL 
POLYNOMIAL  METHOD 

6.65 

0 
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mechods,  with  the  modified  secant  method  selected  three  times  as  often 
as  the  secant  method.  From  these  results,  it  can  be  stated  that  the 
modified  secant  method  introduces  the  lowest  amount  of  error  into  the 
da/dN  data  of  the  six  da/dM  calculation  methods  selected. 

10.3  Distribution  of  da/dN 

The  third  objective  to  be  met  was  to  determine  the  distribution  of 
da/dN  as  a function  of  ^K.  The  first  set  of  da/dN  data  selected  for 
analysis  was  da/dN  data  calculated  by  the  secant  method,  with  the  anti- 
cipation of  also  finding  the  distribution  of  da/dN  data  calculated  by 
the  modified  secant  method  and  the  quadratic  7-point  incremental  poly- 
nomial method.  Data  were  selected  from  the  first  two  da/dN  calculation 
mechods  because  of  their  ability  to  re-create  the  original  a vs.  N data 
and  the  quadratic  7-polnt  incremental  polynomial  method  because  of  its 
widespread  use.  The  combined  data  from  each  of  these  three  methods  are 
shown  in  Figures  36,  37,  and  38. 

The  steps  of  analysis  for  the  distribution  of  da/dN  are  very  similar 
to  the  steps  of  analysis  used  for  Che  distribution  of  N.  First,  the  re- 
plicate da/dN  data  used  was  obtained  from  the  da/dN  vs.  6K.  data  generated 
by  the  DADNCP  program  (Section  7.2)  using  the  secant  method.  Typical  re- 
plicate da/dN  data  are  shown  in  Figure  39.  The  distribution  of  the  re- 
plicate da/dN  data  was  determined  at  each  level  through  the  use  of  the 
CGRDDP  program  (Sectlo  6.3).  At  each  AK  level,  this  program  calculated 
the  distribution  parameters  and  goodness  of  fit  crlr  ~ia  for  the  six  dis- 
tributions and  then  compared  the  goodne  i of  fit  crite.  . between  the  dif- 
ferent distributions  Co  give  Che  distribution  rankings.  Again,  the  gen- 
eralized A-parameCer  gamma  distribution  was  not  included  in  the 
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distribution  rankings.  The  distribution  parameters,  goodness  of  fit 
criteria,  and  the  distribution  rankings  were  then  combined  over  all  of 
the  AK  levels. 

The  distribution  parameCerB  of  the  da/dN  data  as  a function  of  crack 
length  (essentially  ^K)  were  plotted  for  each  of  the  six  distributions 
and  are  shown  in  Figures  40  through  4S.  The  distribution  parameters  are 
again  normalized  to  show  the  trends  present  in  the  parameters. 

The  goodness  of  fit  criteria  for  each  distribution  were  averaged 
over  all  of  the  AK  levels.  Those  results  are  shown  In  Table  IX.  From 
these  results,  an  understanding  of  which  distributions  provide  the  best 
fit  for  the  da/dN  data  can  be  obtained. 

The  distribution  rankings  at  each  AK  level  were  combined  over  all  of 
the  AK  levels  and  again  the  mean  rank  and  Its  standard  deviation  for 
each  of  the  distributions  and  the  number  of  times  each  distribution  was 
selected  as  Che  best  discribuClon  were  calculated.  These  results  are 
shown  in  Table  X. 

Each  of  the  distributions  gives  a fair  but  not  outstanding  performance 
in  providing  a fit  for  the  da/dN  data.  There  were  no  significant  differ- 
ences between  Che  means  of  any  of  the  five  distributions,  especially  con- 
sidering the  high  values  of  standard  deviation  about  the  mean.  The  3- 
parameter  gamna  distribution  did  have  a slightly  lower  mean  than  the 
other  dietribucions  and  it  also  had  the  lowest  value  of  the  Kolmogorov- 
Smlrnov  atatistic.  However,  the  2-peraBister  log  normal  dlatrlbution  was 
the  best  distribution  slightly  more  often  than  the  other  four  distribu- 
tions, but  again  there  were  no  significant  differences  between  the  five 
distributions.  These  results  lead  to  the  conclusion  that  the  3-parainscar 
gamma  distribution  provides  a better  tic  for  the  da/dN  data  than  the  other 
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Figure  43.  3-Paraaieter  Waibull  Diatrlbutlon  Parametare 
of  da/dK  Data  aa  a Function  of  Crack  Length 


NORMALIZED  DISTRIBUTION  PflRDMETERS 


NORMRLIZED  DISTRIBUTION  PARAMETERS 


1.200-1 


l.OOOH 


GENERALIZED  4-PflRAMETER 
GAMMA  DISTRIBUTION 


NORMRLIZED  PRRRMETER 

SECANT  rCTHOD 
DELTA  P = H.20  KIP 
P MAX  = 6.2S  KIP 
DELTA  A = .20  MM 
AO  = 9.00  MM 
NOATA  s 68 

R = .20  X « 


VALUES 

O - TflU  HAT 
A - ALPHA  HAT 
X - B HAT 
^ - G HAT 


Figure  45.  Geoeralitad  4>P«raaMter  Gaaaa  DletribuCion 
Faraneters  of  da/dN  Daca  as  a Function  of 
Crack  Length 


114 


Table  IX 


Average  Goodoeee  of  Fit  Criteria  for  the 
Distribution  of  da/dM  Data 


. CKI-SQUflRE  KOLMOGOROV-  CLOSENESS 

DISTRIBUTION  reiL  arer  smirnov  test  (R  souareoi 


2-PflRRMETER 

NORMAL 

0.8494 

0.0915 

0.94997 

2-PfiRflMETER 
LOG  NORMAL 

0.9011 

0.0779 

0.97647 

3-PARAMETER 
LOG  NORMAL 

0.8442 

0.0834 

0.96966 

3-PflRflMETER 
WEI BULL 

0.8474 

0.0777 

0.95942 

3-PflRflMETER 

GRHMfl 

0.8389 

0.0737 

0.96662 

GENERALIZED 

4-PARAMETER 

GAMMA 

0.7946 

0.0726 
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Table  X 


Distribution  Rankings  for  the  Distribution  of  da/dN  Data 


NUMBER  OF 

STANDARD  TIMES  BEST 

DISTRIBUTION  mean  deviation  distribution 


2-PflRflMETER 

NORMRL 

3.684 

1 .6497 

27 

2-PP)RflMETER 
LOG  NORMAL 

2.603 

1 .1943 

37 

3-PflRflMETFR 
LOG  NORMAL 

3.360 

1 .5524 

26 

3-PARflMETER 
WE I BULL 

2.985 

1.1925 

19 

3-PRRHMETER 

GRMMfl 

2.368 

0.9646 

27 

116 


four  distributions,  but  Its  perforasnca  relative  to  the  other  distribu- 
tions Is  not  strong  at  all.  Due  to  this  poor  performance  by  the  da/dN 
data  in  fitting  a distribution,  no  analysis  of  da/dN  data  calculated  by 
either  the  modified  secant  method  or  the  quadratic  7-polnt  incremental 
polynomial  method  was  conducted. 

10.4  Prediction  of  a vs.  N Data  from  the  Distribution  of  da/dN 

The  fourth  objective  to  be  met  was  to  determine  the  variance  of  a set 
of  a vs.  N data  predicted  from  the  da/dN  distribution  parameters.  The 
da/dN  distribution  parameters  were  estimated  by  the  CGRDDP  program  (Sec- 
tion 6.3)  as  described  in  Section  10.3.  The  AVNFRD  program  (Section  7.3) 
was  run  on  the  da/dN  distribution  parameters  and  68  replicate  a vs.  N 
data  sets  ware  predicted.  These  data  sets  are  shown  in  Figure  46. 

To  obtain  the  variance  of  this  predicted  data,  the  CCDDF  program  (Sec- 
tion 6,2)  was  run  at  14  crack  Length  levels  of  the  predicted  data.  The 
distribution  parameters,  goodness  of  fit  criteria,  and  the  distribution 
rankings  were  then  combined  over  all  of  the  crack  length  levels  run. 

For  this  predicted  data,  neither  the  3-parameter  gaoma  distribution 
nor  the  generalized  4-paraniecer  ganraa  distribution  would  converge  on 
parameter  estimates,  implying  that  neither  distribution  would  provide  a 
fit  for  the  data.  The  distribution  parameters  as  a function  of  crack 
length  obtained  for  the  other  four  distributions  are  shown  in  Figures  47 
through  50.  The  average  goodness  of  fit  criteria  for  the  four  distribu- 
tions for  the  predicted  date  are  shown  in  Table  XI.  The  distribution 
rankings  results  for  the  four  distributions  for  this  data  are  shown  in 
Table  XII. 
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from  the  Distribution  of  da/dN 
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Figure  48.  2'Par«ffl«t«r  Log  Normal  Dlatrlbutlon  Paramatara 
aa  a Funcclon  of  Crack  Langch  for  Cycle  Count 
Data  Pradictad  from  tha  Dlatrlbutlon  of  da/dN 
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Tabl*  XI 

Average  Goodneee  of  Pit  Criteria  for  the  Dlatrlbutlon  of  Cycle 
Count  Data  Predicted  fron  the  Dletrlbutlon  of  da/dM 


CHI-SQ(JfW£  KOLMOSOROV-  CLOSENESS 

DISTRIBUTION  tail  area  snirnov  test  cr  squared) 


2-PflRflMETER 

NORMAL 

0.3087 

0.073S 

0.98515 

2-PflRflMETER 
LOG  NORMAL 

0.9128 

0.0722 

0.98497 

3-PflRflMETER 
LOG  NORMAL 

0.8828 

0.0730 

0.98515 

3-PARRMETER 

WEIBULL 

0.8919 

0.0818 

0.96884 

T«ble  XII 

Dlftrlbution  fUnklogt  for  th«  Distribution  of  Cycle  Count 
Date  Predicted  fron  the  Dietributlon  of  d«/dll 


DISTRIBUTION  mean 


NUMBER  OF 
STRNDf«D  TIMES  BEST 
DEVIATION  DISTRIBUTIt3N 


2-PflRflMETER 

NORMAL 

2.6H3 

0.7449 

2-PflRflMETER 
LOG  NORMAL 

1.214 

0.5789 

3-PflRflMETER 
LOG  NORMAL 

2.286 

0.6112 

3-PARAMETER 

WEIBULL 

3.857 

0.5345 

The  2-paTaiiieter  log  normal  distribution  provided  the  best  fit  for  the 
predicted  replicate  cycle  count  data,  followed  by  the  3>parameter  log 
normal  distribution  and  then  the  2-parameter  normal  distribution.  The  3- 
paraneter  Weibull  distribution  provided  the  worst  fit  for  Che  data  of 
the  four  distributions  which  the  data  fit. 

The  next  step  in  the  analysis  was  the  comparison  of  the  distributions 
of  N between  the  actual  cycle  count  data  and  the  cycle  count  data  pre- 
dicted from  the  distribution  of  da/dN.  The  mean  and  standard  deviation 
of  both  distributions  at  the  crack  length  levels  used  above  were  computed 
and  the  resultc  are  shown  in  Table  aIII.  At  every  crack  length  level, 
there  was  no  significant  difference  between  the  means  but  there  wau  a 
very  significant  difference  between  the  standard  deviations  of  the  two 
distributions.  In  every  case,  the  standard  deviation  of  Che  predicted 
cycle  count  deta  is  much  smaller  than  the  standard  deviation  of  the  actual 
cycle  count  deta. 

As  a check  on  the  analysis  above,  a vs.  N data  were  predicted  from 
the  distribution  of  da/dN  in  a slightly  different  manner  than  for  the 
predicted  replicate  cycle  count  data.  The  mean  and  + 2,  and  3 sigma 

values  of  da/dN  at  each  crack  length  level  were  obtained  from  the  distri- 
bution of  da/dN.  Using  these  7 lines  of  da/dN  data,  a vs.  N data  was 
predicted.  The  results  are  shown  in  Figure  51.  A comparison  between  the 
actual  cycle  count  mean  and  ^ 1,  2,  and  3 sigma  values  and  the  cycle 
count  values  predicted  from  the  mean  and  4 1,  2,  and  3 sigma  da/dN  lines 
at  a single  crack  length  level  is  shown  in  Table  XIV. 

From  the  above  analysis,  it  can  be  concluded  that  predicting  a vs.  N 
data  from  the  distribution  of  da/dN  using  the  method  described  in  Section 
7.2  yields  low  error  in  predicting  mean  crack  propagation  behavior,  but 
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Table  XXII 


CooparLson  of  the  Dletrlbutlons  Betveen  Actual  Cycle  Count  Data  and 
Cycle  Count  Data  Predicted  from  the  Distribution  of  da/dM 


rCAN  STMlOMa  DEUlATIOri 


OMCX  LCWTH  (l«> 

8CTUM. 

nCOXCTCO 

fCTlM. 

PtOlCTED 

11.000 

88681 

88738 

6886 

1991 

13.000 

90186 

91888 

8838 

8918 

IS. 000 

117486 

118700 

6719 

4348 

17.000 

138388 

138871 

10903 

4898 

18.000 

188388 

186918 

11849 

3868 

ai.ooo 

170786 

171379 

18489 

3849 

83.000 

188198 

183670 

8804 

3894 

88.000 

198978 

194804 

8847 

3380 

87.000 

302838 

804083 

8944 

3432 

88.000 

811030 

818618 

9183 

3407 

31.000 

818688 

880383 

9388 

3416 

33.000 

888499 

887186 

9637 

3830 

38.000 

831416 

833849 

10037 

3579 

36.800 

834873 

836533 

10191 

3574 
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T«bl«  XXV 


Conparltotk  of  Xetu«l  Cyol*  Count  Data  vith  Cycle  Count  Data  Predicted 
froa  Conetant  Variance  da/dM  Linea 


- 31. coo  m 

(CTUM. 

neoxcTo 

-3  stem 

18788S 

117430 

-e  stem 

195B8S 

141841 

•1  stena 

aossea 

173198 

NEM1 

ei7991 

813318 

StOM 

^864888 

8S389S 

330341 

♦3  stem 

S7804S 

^ 484830 

yields  high  error  in  predicting  crack  propagation  behavior  at  the  axtreatta 
of  the  distribution  of  N. 

10. S Inverse  Growth  Rata 

Due  to  the  failure  of  Che  da/dH  data  to  fit  any  of  Che  given  distri* 
buttons  satisfactorily,  It  was  decided  that  the  growth  rate  variable  war* 
ranted  a further  investigation.  Looking  back  at  the  original  experloenc- 
al  investigation  (Section  9),  it  can  be  seen  chat  the  dependent  variable 
of  the  data  was  M while  the  independent  variable  was  a (i.e.  N was  mea- 
sured as  a was  varied).  Since  the  dependent  variable,  N,  provided  a vary 
nice  fit  to  the  3-paranieter  log  normal  distribution,  it  was  strongly  sus- 
pected that  changes  in  the  dependent  variable,  ^N,  would  also  provide  a 
good  fit  to  one  of  the  given  distributions.  Since  A a was  constant  It 
was  decided  to  use  AN/da,  or  in  differential  terms,  dN/da,  as  a variable 
of  interest  for  further  analysis. 

The  analysis  conducted  using  dN/da  as  the  variable  of  interest  was 
Che  same  analysis  used  for  da/dN.  The  first  part  of  this  analysis  was  to 
determine  the  distribution  of  dN/da.  Replicate  dN/da  data  were  obtained 
by  inverting  Che  replicate  da/dN  dsta  calculated  by  the  secant  method 
using  Che  DADNCF  program  (Section  7.2).  Typical  replicate  dN/ds  data  are 
shown  in  Figure  52.  The  distribution  of  Che  replicate  dN/da  data  was 
determined  at  each  AK  level  through  the  use  of  the  DMDDF  program  (Section 
6.4).  At  each  AK  level,  this  program  calculated  the  distribution  para- 
maters  and  goodness  of  fit  criteria  for  six  distributions  and  than  com- 
pared Che  goodness  of  fit  criteria  bpcwean  Che  different  distributions 
to  give  the  distribution  rankings.  The  location  parameter  for  the  3- 
parameter  gamma  distribution  and  the  generalised  4-parsmster  gamma 
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dlscrlbutlon  was  assumed  co  be  eero  by  this  program,  thus  reducing  these 
two  distributions  to  the  2*paraneter  gamma  distribution  and  the  general- 
ised 3-paraineter  gamma  distribution,  respectively.  As  before,  the  gen- 
eralised 3-parameter  gamma  distribution  was  not  included  in  the  distribu- 
tion rankings.  The  distribution  parameters,  goodnesa  of  fit  criteria, 
and  the  distribution  rankings  were  then  combined  over  all  of  the  AK  levels. 

The  distribution  parameters  of  the  dN/da  data  as  s function  of  crack 
length  were  plotted  for  each  of  the  six  distributions  and  are  shown  in 
Figures  53  through  58.  The  distributiorfs  parameters  are  again  normalised 
to  show  the  trends  present  in  the  parameters. 

The  goodness  of  fit  criteria  for  each  distribution  were  averaged  over 
all  of  the  AK  levels.  These  results  are  shown  in  Table  XV.  From  these 
results,  an  understanding  of  which  distributions  provide  the  best  fit  for 
the  dN/da  data  can  be  obtained. 

The  distribution  rankings  at  each  AK  level  were  combined  over  all  of 
the  AK  levels  and  the  mean  rank  and  its  standard  deviation  for  each  of 
the  distributions  and  the  number  of  times  each  distribution  was  selected 
as  the  best  distribution  were  calculated.  These  results  are  shown  in 
Table  XVI. 

The  3-paraffleter  log  normal  distribution  provided  the  best  fit  for  the 
dN/da  data,  as  evidenced  by  the  low  distribution  ranking,  the  low 
Kolmogorov-Smirnov  test  statistic  value,  and  the  large  number  of  times 
it  was  selected  as  the  best  distribution.  The  2-paraBieter  log  normal  and 
the  3-parameter  Weibull  dietribution  tied  for  the  second  best  fit  for  the 
dM/da  data,  both  having  roughly  the  same  distribution  ranking  and 
Kolmogorov-Smirnov  test  statistic  value  and  the  same  number  of  timos  it  was 
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Figure  bA.  2*Paraaieter  Log  Normal  Discributlon  ParaoMtara 
of  dN/da  Data  aa  a FuncClon  of  Crack  length 
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TAbU  XVI 

Distribution  KAnklngs  for  ch«  Distribution  of  dN/da  Data 
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vslccMd  a>  tha  bast  dlatrlbutloa.  The  2-paranaCar  ganma  distribution 
provided  tha  fourth  bast  fit  aod  tha  S^paraiaetar  normal  distribution  pro- 
vided the  worst  fit  for  tha  dN/da  data. 

The  next  step  of  tha  analysis  was  to  see  if  the  Ingiroved  fit  of  the 
dJi^da  dafea>  to  a distribution  would  improve  tha  prediction  of  a vs,  N data 
from  the  distribution  of  dN/da.  The  AVimD  program  (Section  7.3)  was 
slightly  modified  for  tha  dN/da  variable  and  run  on  the  dN/da  distribution 
paramatarS(  resulting  in  68  predicted  replicate  a vs.  N data  sets.  These 
data  seta  are  shown  in  Figure  59. 

The  CCDDP  program  (Section  6.2)  was  run  at  a few  crack  length  levels 
of  the  predicted  data.  The  distribution  parameters,  goodness  of  fit 
criteria,  and  the  distribution  ranki.igs  were  then  combined  over  all  of  the 
crack  length  levels  run. 

For  this  set  of  predicted  data,  the  generalised  4-paramBter  games 
distribution  would  not  converge  on  parameter  estimates,  implying  that  it 
could  net  provide  a fit  for  the  dN/de  data.  The  distribution  paramtors 
as  s function  of  crack  length  obtained  for  the  other  five  distrlbutiona 
are  shown  in  Figures  60  through  64.  The  average  goodness  of  fit  critsria 
for  ths  five  distributions  for  the  predicted  data  ere  ehown  in  Table  XVII. 
The  distribution  rankings  results  for  the  five  distributions  fpr  this  data 
are  shown  in  Table  XVIII, 

The  3-paramstsr  log  normal  distribution  provided  the  best  fit  for  tha 
pradictsd  repllcsCa  cycle-count  data,  followed  in  order  by  the  2-paraBwtaf 
log  normal  distribution,  ths  3-paraineter  Weibull  distribution,  ths  2-psrS' 
meter  normal  distribution,  sod  the  3-parsaeter  gamma  distribution. 

Ths  next  step  in  tha  analysis  was  the  comparison  of  the  distributions 
of  N between  the  actual  cycle  count  data  aod  the  cycle  count  data  predicted 
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Table  XVII 


Average  Goodness  of  Fit  Criteria  for  the  Distribution  of  Cycle 
Count  Data  Predicted  from  the  Distribution  of  dN/da 


CHI -SQUARE  KOLMOGOROV-  CLOSENESS 
DISTRIBUTION  TAIL  area  SMIRNOV  test  IR  SQUARED) 


2-PflRflMETER 

NORMAL 
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0.0652 

0.98765 
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LOG  NORMRL 
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LOG  NORMAL 
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Table  XVIIl 


Distribution  Rankings  for  the  Distribution  of  Cycle  Count  Data 
Predicted  from  the  Distribution  of  dN/da 


DISTRIBUTION  MEAN 


NUMBER  OF 
STANDARD  TIMES  BEST 
DEVIATION  DISTRIBUT™ 


from  the  distribution  of  dN/da.  The  mean  and  standard  deviation  of  both 
distributions  at  the  crack  length  levels  used  above  were  computed  and  the 
results  are  shown  in  Table  XIX.  At  every  crack  length  level,  there  was 
no  significant  difference  between  the  means  but  there  was  a very  signifi- 
cant difference  between  Che  standard  deviations  of  the  two  distributions. 

In  every  case,  the  standard  deviation  of  the  predicted  cycle  count  data 
was  much  smaller  thai.  the  standard  deviation  of  the  actual  cycle  count 
data. 

Just  as  for  Che  data  predicted  from  the  distribution  of  da/dN,  a vs. 

N data  were  predicted  from  Che  mean  and  + 1,  2,  and  3 sigma  dN/da  lines. 

The  results  are  shown  in  Figure  65.  A comparison  between  the  actual  cycle 
count  mean  and  1,  2,  and  3 sigma  values  and  the  cycle  count  values  pre- 
dicted from  Che  mean  and  +1,  2,  and  3 sigma  dN/da  lines  at  a single  crack 
length  level  is  shown  in  Table  XX. 

From  the  above  analysis,  it  can  be  concluded  again  that  predicting  a 
vs.  N data  from  the  distribution  of  dN/da  using  Che  method  described  in 
Section  7.3  yields  low  error  in  predicting  mean  crack  propagation  behavior, 
but  yields  high  error  in  predicting  crack  (propagation  behavior  at  the  ex- 
tremes of  the  distribution  of  N. 


Table  XIX 


I a 

• ] 


Comparison  of  Che  Dlscrlbutlona  Between  Actual  Cycle  Count  Data  and  ! ^ 

Cycle  Count  Data  Predicted  from  the  Distribution  of  dN/da  i ^ 
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Comparison  of  Actual  Cycle  Count  Data  with  Cycle  Count  Data 
Predicted  from  Constant  Variance  dN/da  Lines 
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DISCUSSION 

Throughout  the  course  of  this  investigation,  a few  unique  events  tojk 
place  and  many  interesting  observations  were  made.  Some  of  these  have 
rather  simple  explanations  r \ others  require  a quite  detailed  discus- 
sion. Hopefully,  some  important  conclusions  can  be  made  as  a result. 

11,1  Experimental  Investigation 

The  behavior  of  fatigue  crack  propagation  experienced  during  this 
investigation  was  much  different  than  first  anticipated.  The  most  sur- 
prising event  that  took  place  in  almost  every  test  was  the  sudden  changes 
in  the  magnitude  of  the  crack  growth  rate.  Both  sudden  increases  in  the 
growth  rate,  as  if  the  crack  had  just  cooae  upon  some  unusually  weak  aIu^^ 
inum,  and  sudden  decreases  in  the  growth  rate,  as  if  the  crack  was  exper- 
iencing some  unusually  cough  material,  were  observed  repeatedly,  many 
times  one  or  two  millimeters  after  a previous  event  of  similar  nature. 

One  of  the  more  outstanding  examples  cf  this  type  of  behavior  is  shown  in 
Figure  66. 

It  appears  chat  the  material  is  made  up  mostly  of  a fairly  homogeneous 
material  with  many  smaller  areas  located  in  a random  fashion  which  char- 
acteristically have  vastly  different  crack  propagation  properties  Chan  Che 
majority  of  the  material.  The  size  of  these  small  areas  seems  to  vary 
considerably  from  as  small  as  less  Chan  1 millimeter  In  length  to  perhaps 
as  large  as  S or  10  mlllimetera  in  length.  These  small  areas  obviously 
have  a very  large  effect  on  Che  overall  smoothness  of  an  a vs.  N data  set 
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and  on  the  total  amount  of  scatter,  especially  in  the  growth  rate  data. 

The  average  growth  rate  also  seemed  to  vary  somewhat  from  test  Co 
test,  with  some  tests  running  slow  throu^out  Che  whole  test,  while 
other  tests  ran  fairly  fast  throughout  the  whole  test.  This  phenomenon 
la  Che  cause  of  Che  outlying  data  sets  in  Figure  22.  As  also  noted  by 
ocher  Investigators  [5,6,7],  the  variation  in  growth  rate  at  the  beginning 
of  the  test  during  Che  slow  growth  rates  leads  to  most  of  the  variation 
In  N at  the  final  crack  length. 

As  a result  of  these  observations,  the  conclusion  Is  made  that  this  al- 
loy Is  a very  non-homogeneous  material,  especially  considering  the  random 
nature  of  crack  propagation  behavior.  It  very  rarely  obeys  the  smooth 
growth  rats  equations  often  used  to  describe  Its  behavior  and  does  so 
only  when  it's  behavior  Is  considered  at  a very  macroscopic  level. 

11.2  Distribution  of  N 

The  conclusion  stating  chat  the  cycle  count  data  follows  a 3-paraniecer 
log  normal  distribution  can  be  considered  very  strong.  The  only  occur- 
rences where  chjs  was  not  so  was  at  short  crack  lengths  where  the  need 
for  the  location  parameter  used  in  the  S-parameter  log  normal  distribu- 
tion was  not  near  as  strong  as  at  long  crack  lengths. 

Fromflgures  24  through  29,  the  distribution  parameters  tend  to  vary 
quite  a bit  at  short  crack  lengths  but  ten',  to  follow  smooth  curves  after 
a ^ IS  mm.  The  scale  parameter  In  the  first  two  distributions  where  no 
location  parameter  Is  estimated  have  very  smooth  curves,  showing  that 
mean  crack  propagation  behavior  does  follow  smooth  growth  rate  equations. 
The  sane  smooth  shape  of  the  location  pramster  In  ths  last  four  discrlbu- 
’tlons  also  supports  this  Staten. nt.  Essentially  these  location  paraneter 
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curves  define  aw  area  where  crack  propagation  will  never  occur.  In  other 
words,  the  number  of  load  cycles  needed  to  reach  a given  crack  length 
will  never  be  less  than  the  estimate  of  the  location  parameter  at  that 
crack  length.  This  is  shown  in  Figure  67.  On  this  plot,  crack  propaga- 
tion data  will  never  occur  to  the  left  of  the  location  parameter  line. 

Note  from  Figure  26  that  the  scale  parameter,  tends  to  remain  constant 
after  a ^ 20  am.,  allowing  the  location  parameter  to  completely  account 
for  the  increase  of  H with  increasing  crack  length.  From  Figure  24,  note 
the  smooth  increase  of  the  shape  parameter,  a , as  a function  of  crack 
length,  thus  supporting  the  expectation  of  higher  variances  at  longer 
crack  lengths.  Another  interesting  event  is  shown  in  Figure  29.  The 
power  parameter,  o,  of  the  generalized  4-porameter  gamma  distribution  was 
always  estimated  to  be  equal  to  one,  thus  reducing  this  distribution  to 
the  3-paramecer  gamma  distribution.  It  should  be  noted  here  that  over 
half  of  the  computer  time  used  to  obtain  all  of  the  distribution  para- 
meters was  used  to  estimate  the  parameters  of  the  generalized  4-parameter 
gamma  distribution.  By  eliminating  this  distribution  from  the  CCDDP  pro- 
gram, much  time  and  money  can  be  saved. 

Another  interesting  occurrence  which  appeared  very  often  is  shown  in 
Table  VI.  Many  times  the  distribution  rankings  implied  by  one  goodness 
of  fit  criterion  could  not  be  supported  by  another  goodness  of  fit  cri- 
terion and  often  three  different  distribution  rankings  were  implied  by  the 
three  goodness  of  fit  criteria.  In  other  words,  the  goodness  of  fit  cri- 
teria were  not  consistent  between  themselves.  This  necessitated  a some- 
what subjective  analysis  of  the  goodness  of  fit  criteria.  The  closeness, 

2 

R , tended  to  be  very  sensitive  to  the  scales  of  the  plot  and  the  slope 
of  the  linear  least  squares  line.  Thus,  the  closeness  was  rarely  used 
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Figure  67.  Eecloiace  of  che  Location  Parameter  of  the 
3>Farameter  Log  Normal  Diatrlbutlon  aa  a 
Function  of  Crack  Length 


unless  the  slopes  were  approximately  the  same  between  the  different  dis- 
tributions. The  chl-squarc  tail  area  tended  to  be  undlscernlng  between 
distributions  that  provided  fairly  equal  fits  to  the  data.  Thus,  the  chi- 
square  goodness  of  fit  criterion  vae  used  only  when  there  were  fairly 
large  differences  between  the  distributions.  The  Rolmogorov-Smlrnov 
goodness  of  fit  teat  provided  a fairly  reliable  and  sensitive  test  and 
was  used  heavily  in  establishing  distribution  rankings. 

A typical  fit  provided  by  each  of  the  five  distributions  for  the 
cycle  count  data  is  shown  in  Figures  68  through  72.  As  state  previously, 
the  3-parameter  log  normal  distribution  provided  a reliable  tight  fit  for 
the  cycle  count  data  as  shown  In  Figure  70.  The  3- parameter  gamma  dis- 
tribution did  surprisingly  well  and  although  it  was  not  selected  as  the 
best  distribution  very  often,  it  consistently  placed  a close  second  to 
the  3-paraoeter  log  normal  distribution.  The  2-paraiiieter  log  normal  dis- 
tribution did  not  do  well  due  to  the  lack  of  a location  parameter.  For 
the  3-parameter  Weibull  distribution,  the  location  parameter  seemed  to 
work  alright,  but  the  shape  of  the  density  function  did  not  match  the  data 
very  well  as  shown  in  Figure  71.  The  2-parameter  normal  distribution 
provided  a very  poor  fit  for  the  cycle  count  data  and  should  not  be  in- 
cluded in  any  further  investigations  of  the  dictribution  of  N. 

11.3  Crack  Growth  Rate  Calculation  Methods 

Of  the  six  da/dN  calculation  methods  selected,  both  the  secant  method 
and  the  modified  secant  method  contributed  low  amounts  of  error  into  the 
da/dN  data  as  shown  in  Table  VIII.  The  modified  secant  method  calculated 
da/dN  data  which  could  be  Integrated  back  closer  to  the  original  a vs.  N 
data  chan  the  secant  method  could,  perhaps  because  it  calculates  da/dN 
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Figure  68.  Typlsal  Fic  of  the  Cycle  Count  Date  to 
the  2'Paraaeter  Normal  Dlatrlbutlon 
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Figure  70.  lypicel  Fit  of  the  Cycle  Count  Date  to 
the  3*P«reffleter  Log  Norael  Dietrlbutlon 
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Figure  72.  Typical  Fit  of  the  Cycle  Count  Data  to 
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d^ta  ar  the  original  crack  length  levels  Instead  of  between  Che  original 
crack  length  levels. 

None  of  th^  InQremeutal  polynomial  ooChods  ca).culaced  da/dN  data 
which  could  be  integ..aced  bc%k  even  close  to  the  original  a vs.  N data. 

This  it.  no  doubt  due  to  the  soo«.>>ing  effect  of  these  methods  which  tends 
to  reduce  the  sudden  changes  in  growth  rates.  This  is  shown  in  Figure  38 
where  the  number  of  extreme  da/dN  data  points  tor  the  quadratic  7-point 
incremental  polynomial  method  is  much  less  than  the  number  of  extreiM 
da/dN  data  points  for  either  Che  secant  method  of  Che  modified  secant 
method  (Figures  36  and  37).  This  is  also  shown  in  Figures  30  through  35 
where  the  incremental  polynomial  method  data  lullow  a narrow  band  .line 
while  Che  recant  method  and  modified  secant  method  data  follow  a more 
broad  band  line.  Note  also  from  Figures  32  through  35  chj  waviness  of 
Che  data  show-i.p.g  the  Large  chl^ges  in  growth  rate  noticed  during  data  ac- 
quisition. . . 

■'1 

If  crack  propagation  data  w^e  always  very  smooth  data,  chan  cha  io- 
cramental  polynomial  methods  would  introduce  a very  small  amount  of  error 
into  Che  da/dN  data.  But  as  stated  pt^viously,  the  sudden  changes  in 
growth  rata  are  inherent  in  the  crack  propagation  process,  and  any  attaopt 
at  modifying  these  changes  will  distort  tbf  resulting  date  end  prevent  It 
from  becoming  e true  representetion  of  crack  'pafopegati:>n  behavior. 

Of  the  four  Incremental  polynomial  methods  used,  both  bbe  gu^ratlc 
7-polnt  version  and  the  qu^retlo  log-leg  7-fleiat  ^uerslbiil  db  the  best  job, 
followed  closely  by  the  linear  7-polnt  version.  The  linear  log- log  7- 
polnt  version  does  a vary  poor  job  ps  showp  In  Figure  34  and  Table  Vlll. 

The  use  of  the  log- log  transformitlon  failed  to  give  any  Improvamsnt  over 
the  conventional  Incremental  polynomial  methods  In  the  ability  to  reproduce 
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the  originel  a vs.  N data.  The  use  of  the  second  order  polynomial  fit 
over  a straight  line  polynomial  fit  improves  the  performance  of  Che  in* 
cremental  polynomial  meth^s,  especially  when  using  the  log-log  transfor- 
mation. - 

The  amount  of  variaeibn  -ih  the  fverggie  tncreiiien]||[!l  .j^r  cent  error  over 
all  of  tha'^esth  mae  fairly  Iti^ll,  usualiy  leas  than^b  per  cent  error, 
indicating  fairly  consistent  results  over  all  of  Che  experimental  tests. 


11.4  Distribution  of  da/dN 

No  outstanding  positive  results  were  achieved  for  the  diatlrflfittion  of 

’ i 

da/dN.  Each  of  the  distributions  provided  roughly  the  same  quejlity  of  fit 

I 

for  the  data,  with  the  3-parameCer  gaiuna  distribution  doing  a ^lightly 
better  Job  than  the  other  four  distributions.  r 

J 

The  da/dN  data  varied  quite  a bJLc  as,  a function  of  ilK  as  ijhown  in 
Figure  36.  As  a result,  different  distributions  would  provide  |a  fit  for 
Che  da/dN  data  at  different  dK  levels,  depending  on  the  generalT^fS^  and 
skewness  of  the  data  at  a given  dK  level.  There  were  several  (^cas^ons 
when  Che  da/dN  data  was  skewed  leftf,  as  shown  in  Figure  73,  synmcric,  as 
shown  in  Figure  74,  or  skewed  right,  as  shown  in  Figure  75. 

When  Che  data  was  either  skewed  left  or  symmetric,  the  2'iMcameCer 

. i 

normal  distribution  provided  the  beat  fit  for  the  da/dN  data,  ais  shown  in 


Figures  76  and  77.  When  th^  data  waa.ajcpwed  right,  either  log*  nor- 

mal distributions,  the 'S^paramster  Welbbll  dtsrfl'butlon,  or  the  3-parameter 
gansu  distribution  provided  a fit  for  Che  da/dN  data.  Typical  fits  of  the 
skewed  right  da/dN  (fata  tb  thkte  fbur  disttlbutious  iA’l''4(ii&wn  in  Figures 
78  through  81.  Due  to  the  large  variations  in  the  da/dN  data,  each  of  the 
distributions  is  needed  to  provide  a fit  for  the  wide  range  of  density 
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function  shapes. 

The  large  amount  of  variation  in  the  shape  of  the  da/dN  data  aa  a 
function  of  AK  is  shown  In  the  plots  of  the  distribution  parameters  as  a 
function  of  crack  length  (Figures  40  through  45).  Note  that  there  Is  a 
lot  of  variation  indie  shape  parameter  (the  pluses)  and  the  location  para> 
meter  (the  diamonds).  The  variation  of  the  shape  parameter  reflects  the 
changes  in  the  amount  of  variance  and  the  shape  of  the  data.  The  variation 
of  the  location  paraiseter  reflects  the  changing  skewness  of  the  data.  As 
the  skewness  goes  from  right  to  left,  the  estimate  of  the  location  para* 
meter  decreases  rapidly.  Also*  from  Figure  45,  it  can  be  seen  that  there 
are  many  occurrences  where  the  estimate  of  a was  not  equal  to  one,  thus 
Implying  the  necessity  of  the  inclusion  of  the  generalized  4-parameter 
game  distribution  when  analysing  da/dN  data  so  that  a wide  range  oi  den- 
sity function  shapes  can  be  acconuodated  for  the  da/dN  data. 

11.5  Prediction  of  a vs,  N Data  from  the  Distribution  of  da/dN 

The  results  of  the  prediction  of  replicate  a vs.  N data  from  the  dis- 
tribution of  da/dN  were  less  revealing  than  anticipated.  When  comparing 
Figure  46  with  Figure  22,  it  becomes  apparent  that  the  variance  of  the 
predicted  a vs.  M data  Is  much  less  than  the  variance  of  the  actual  a vs. 

N data.  However,  the  mean  of  the  predicted  a vs.  N data  is  very  close  to 
the  mean  of  the  actual  a vs.  N data.  The  implication  of  this  is  that 
crack  propagation  behavior  is  not  being  accurately  modeled  by  a randomly 
selected  value  of  da/dN  from  the  distribution  of  da/dN,  In  crack  propa- 
gation behavior,  as  discussed  in  Section  11.1,  the  growth  rate  at  a given 
AK  level  Is  not  independent  of  the  growth  rates  at  previous  AK  levels,  as 
evidenced  by  periods  of  up  to  10  ism.  of  uncharacteristically  fast  or  slow 
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growth  rates.  However,  the  independence  of  growth  rates  is  assumed  in 
the  prediction  of  a vs.  N data  from  the  da/dN  distribution  parameters, 
resulting  in  very  sisooth  a vs.  N data.  This  ssmoth  a vs.  N data  lacks 
the  areas  of  sudden  fast  and  slow  growth  rates  discussed  in  Section  11.1 
which  occurs  frequently  in  actual  a vs.  N dcta.  Thus,  the  co«d>ination  of 
isany  smooth  a vs.  H lines  of  the  same  mean  behavior  results  In  Hte  reduc- 
tion of  variance  noted  above.  To  accurately  predict  crack  propagation 
behavior,  some  means  of  quantitatively  describing  the  interdependence  of 
adjacent  growth  rates  must  be  found. 

When  the  distribution  of  the  cycle  count  data  predicted  from  the 
distribution  of  da/dN  was  analyzed,  neither  gamma  distribution  would 
converge  on  its  parameters  as  the  estimate  of  the  shape/power  parameter, 
g,  tended  to  approach  its  upper  global  limit.  The  2-paraineter  log  nonuil 
distribution  provided  the  best  fit  for  this  data  because  the  location 
parameter  of  the  3-parameter  log  normal  distribution  was  estimated  to  be 
zero  at  most  crack  length  levels. 

When  the  distribution  of  the  predicted  cycle  count  data  is  compared 
with  Che  distribution  of  Che  actual  cycle  count  data,  ae  shown  in  Table 
XIII,  it  can  again  be  seen  how  the  mean  of  the  predicted  cycle  count  data 
is  very  close  to  the  mean  of  Che  actual  cycle  count  data  while  the  stan- 
dard deviation  of  the  predicted  cycle  count  data  is  much  less  than  the 
standard  deviation  of  the  actual  cycle  count  data. 

When  a vs.  N data  are  predicted  from  constant  variance  da/dN  lines, 
the  spread  of  the  predicted  data  is  siuch  wider  than  the  spread  of  the 
actual  data,  as  shown  in  Table  XXV.  This  occurs  because  either  all  very 
slow  or  very  fast  growth  rate  data  is  used  at  the  ± 3 sigma  da/dN  lines, 
thus  causing  either  a very  long  or  very  short  number  of  cycles.  The 


176 


actual  data,  however,  rarely  has  any  growth  rates  on  the  order  of  + 3 
Blgna,  and  even  more  rarely  has  repeated  growth  rates  on  the  order  of  + 3 
sigma.  On  the  average,  actual  data  tend  to  have  repeated  growth  rates 
within  jh  1 sigota. 

From  Figure  51,  It  can  be  seen  that  the  constant  variance  lines  tend 
CO  gat  further  apart  when  going  from  left  to  right.  Indicating  that  the 
distribution  of  N is  skewed  right.  Since  the  distribution  of  N has  been 
determined  as  the  3-paraineter  log  normal  distribution  which  is  a skewed 
right  distribution,  the  prediction  of  a vs.  N data  from  constant  variance 
da/dN  lines  supports  this  conclusion. 

11.6  Inverse  Growth  Rate 

As  anticipated,  an  improvement  in  the  fit  provided  for  Che  dN/da 
data  over  the  fit  provided  for  the  da/dN  data  was  obtained.  The  3>para> 
mater  log  normal  distribution  was  able  to  provide  the  best  fit  for  the 
dN/da  data  without  serious  competition  from  the  other  four  distributions. 
This  isiprovement  Is  partially  due  to  the  Inversion  of  the  growth  rate 
variable.  Since  N was  strongly  log  normally  distributed,  it  was  antici- 
pated that  AN  would  be  log  normally  distributed  also.  Another  reason  for 
this  improvemsnt  was  Che  exclusion  of  the  location  paraswCer  from  the 
gaama  distributions,  thereby  severly  decreasing  their  ability  to  provide 
an  adequate  fit  for  the  dN/da  data.  The  fit  provided  by  those  distribu- 
tions which  estimated  a location  paramstar  was  significantly  better  chan 
the  fit  provided  by  the  gasna  distributions.  Quite  a large  range  of 
values  were  estimated  for  the  Location  paramstar  (from  - 1.6  x 10*^^  to 
4.9  X 10^)  and  the  absolute  value  of  the  astiaate  of  tba  location  para- 
■scar  was  always  greater  than  900,  iadicating  no  tendency  to  approach 
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caro  as  acsuned  In  equation  48  (Section  5.2.e)<  Thus,  this  assumption 
has  not  proven  valid. 

The  value  of  the  location  parameter  of  the  3-paraaieter  log  normal  dis* 
trlbution  assumed  negative  values  in  many  instances,  indicating  that 
skewed  left  and  symostrlc  dN/da  data  was  present  as  well  as  skewed  right 
dN/da  data.  This  was  expected  since  the  simple  inversion  of  the  da/dN 
variable  does  nothing  to  change  the  skewness  of  the  density  function  of 
the  data.  The  only  effect  of  this  Inversion  is  to  change  the  direction 
of  the  skewness  and  to  alter  the  shape  of  the  density  function  slightly. 

A histogram  of  typical  symmetric  dN/da  data  is  shown  in  Figure  82  and 
plots  of  the  fit  of  Che  dN/da  data  to  each  of  Che  distributions  are  shown 
in  Figures  83  through  87.  Note  in  Figure  85  Che  ability  of  the  3*para> 
meter  log  norgtal  distribution  to  handle  symnetric  as  well  as  skewed  right 
data.  Again,  due  to  the  large  variation  in  the  dN/da  data,  each  of  the 
distributions  is  needed  in  order  Co  provide  a fit  for  Che  data. 

There  is  a large  variation  in  the  shape  parameter  and  location  para- 
meter again  for  the  dN/da  data,  as  shown  in  the  plots  of  the  distribution 
parameters  as  a function  of  crack  length  (Figures  53  through  58).  The  use 
of  dN/da  does  not  remove  these  variations  from  the  data,  although  it  does 
reverse  the  basic  trend  of  the  mean  as  shown  by  comparing  Figure  40  with 
Figure  53.  The  Mao  value  of  ds/dN  increases  as  a function  of  crack 
length  while  the  nan  value  of  dN/da  decreases  as  a function  of  crack 
length,  both  being  expected  for  constant  amplitude  loading. 

The  use  of  dN/da  distribution  parameters  in  the  prediction  of  repli- 
cate a vs.  N data  did  not  change  the  predicted  data  notloeably.  As  sug- 
gested fxrevlously,  the  problem  of  predicting  a vs.  N data  accurately  lies 


178 


RELATIVE  FREQUENCY 


.2400-1 


.2000 H 


RELfiTIVE 

FREQUENCY  HISTOGRAM 

REPLICATE  CA  TESTS. 

SECANT  HETHOD 

DN/DA  CLASS  SIZE  = 106H4  R = .20 

DELTA  P = H.20  KIP  A = 21.10  MM 

68  DATA  POINTS  DELTA  A = .20  MM 


.leooH 


.1200  H 


.0800 H 


I 


ICYCi^/IN)  •) 


Figure  62.  typical  By— aerie  dH/da  Data 


CUMULATIVE  STONOflRO  NORMAL  PROBABILITY  (PER  CENT) 


2-PRRflMETER 

NORMRL  DISTRIBUTION  PLOT 


DEVIATION  OF  DN/DA  FROM  THE  MEAN 


Flfur*  83.  Typieal  Fit  of  SjMMtrle  81i/d«  Ooto  to  tbo 
2«PttraMt«r  IoxmI  Dlttrlbutioa 


180 


CUMULATIVE  STANDARD  NORMAL  PROBABILITY  (PER  CENT) 
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181 


CUMULfiTIVE  STANDARD  NORMAL  PROBABILITY  (PER  CENT) 
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not  in  which  variable  is  used  to  predict  the  data  but  rather  in  the 
assumption  of  independent  adjacent  growth  rates. 

When  the  cycle  count  data  was  predicted  from  the  dN/da  distribution 
parameters,  the  3>parameter  log  normal  distribution  provided  the  best  fit 
for  the  data  as  the  estimatea  of  the  location  parameter  were  all  at  anti- 
cipated values.  This  is  an  improvement  over  the  prediction  of  cycle 
count  data  from  the  distribution  parameters  of  da/dN,  because  the 
estimate  of  the  location  parameter  was  nearly  always  equal  to  zero.  The 
use  of  this  location  parameter  significantly  Improves  the  fit  of  the  pre- 
dicted cycle  count  data  to  the  3-paraoeter  log  normal  distribution. 

Again,  the  values  of  g assumed  maximum  global  values  in  both  gamma  distri- 
butions. This  is  most  likely  due  to  a lack  of  significant  variance  in 
Che  predicted  cycle  count  data.  When  the  distribution  of  predicted  cycle 
count  data  was  compared  again  to  the  distribution  of  actual  cycle  count 
data,  the  mean  data  was  almost  exactly  predicted  while  the  predicted 
standard  deviation  was  again  much  less  chan  the  actual  standard  deviation, 
which  can  be  seen  by  comparing  Figure  39  with  Figure  22. 

The  a vs.  N data  predicted  from  constant  variance  dM/da  lines  almost 
exactly  reproduced  a similar  plot  made  from  constant  variance  da/dN  lines, 
as  seen  by  comparing  Figure  31  with  Figure  65.  Thus,  the  dN/da  data  seems 
to  support  the  conclusion  that  the  cycle  count  data  fits  the  3-parameter 
log  normal  distribution  the  best. 
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SECTXOK  XII 


CONCLItSlONS 


The  most  significant  conclusions  of  this  Investigation  are  sumnarited 
as  follows: 

1)  The  2>paraiiieter  Welbull  distribution  was  tried  on  previously 
generated  fatigue  crack  propagation  data  and,  due  to  lea  very 
poor  perf oriaance , was  dropped  from  the  remainder  of  the  statisti- 
cal analysis  (Section  8.1) 

2)  Actual  replicate  cycle  count  data  followed  a 3-parameter  log  nor- 
mal distribution,  with  especially  good  fits  at  longer  crack 
lengths  (Section  11.2). 

3)  The  modified  secant  method  Introduces  the  lowest  amount  of  error 
Into  the  da/dN  data  of  the  six  growth  rate  calculation  methods 


selected  (Section  11.3). 

4)  The  large  amount  of  variance  present  In  the  da/dN  vs.  dK  data 
prevented  a consistent  fit  of  the  replicate  da/dN  data  to  any  of 
the  candidate  distributions  (Section  11.4). 

5)  Replicate  dN/da  data  followed  a 3-paramater  log  normal  distribu- 
tion (Section  11.6). 

6)  The  method  of  predicting  a vs.  N data  from  the  da/dN  or  dN/da 


distribution  parameters  was  not  completely  successful  due  to  the 
assumption  of  Independent  adjacent  growth  rates  (Sections  11. S 
and  11.6). 
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RBCOMCMDATlOltS  FOR  FUltTHSR  WORK 

The  use  of  etatlstical  mechoda  In  deacrlMng  aad  predicting  fatlgxie 
crack  propagation  behavior  worked  very  well.  However,  accurate  life 
prediction  waa  not  achieved  becauae  a total  etatlatlcal  deacrlptlon  of 
the  crack  propagation  process  has  not  been  determined.  Only  a minute 
percentage  of  the  total  possible  experimental  and  analytical  work  needed 
to  achieve  this  total  statistical  description  was  conducted  under  this 
investigation.  Based  on  the  observations,  results,  and  conclusions  of 
this  investigation,  the  following  topics  need  further  Investigation. 

1)  Experimental  crack  propagation  data  with  M as  the  independent 
variable  and  a as  the  dependent  variable  Is  needed.  From  this 
the  distribution  of  a as  a function  of  N and  the  distribution  of 
da/dN  as  a function  of  N can  be  obtained. 

2)  A study  of  the  Interdependence  of  growth  rate  data  would  be 
valuable  for  use  In  the  prediction  of  a vs.  N data  from  the  dis- 
tribution of  growth  rate  data. 

3)  The  affect  of  data  density  on  the  variance  and  distribution  of 
growth  rate  data  needs  to  be  found  to  eld  In  more  accurate  data 
acquisition  and  analysis. 

4)  A study  of  the  sudden  growth  rate  changes  In  the  original  a vs.  N 
data  mantlonsd  In  Section  11.1  would  aid  considerably  In  the 
understanding  of  the  crack  propagation  process. 
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S)  A more  reliable  and  accurate  method  of  eatabllshlng  the  distri- 
bution rankings  la  needed.  The  goodness  of  fit  criteria  used  in 
this  investigation  did  not  totally  fulfill  this  ,naad. 
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APPENDIX  A 


D6RIVATZ0N  OP  THE  da/dN  EQUATION  FOR  THE  LZIBAR  LOG-ZX)G  7-POZNE 
INCRBWNIAL  POLYNOMIAL  METHOD 


The  fitted  polyoomlal  equation  for  the  linear  log-log  7-point  Incre- 
oental  polynomial  method  ie  given  by 


'>1  "is 


(A-1) 


where  N..  ie  given  by 
L9 


^ - "l 


(A-2) 


where  and  ^re  given  by  the  ecaling  aquatione  (aquationc  7 and  8, 
Section  4.3).  Subetituting  into  equation  A-1  for  *4.8» 


r^®*10  ^ 

‘“»io  • • ‘o  * ‘i  L — c; — 


Solving  for  a, 


a ■ 10 


b + b, 
o 1 


N - C, 


b,  C,  b,  log,„  N 


’’o  Cj 


"V  ““‘lO  ">  • cj 
^ . 10  ^ 


**0  ■ C 


b,  C,  b. 


o Cj  Cj 
- 10  * • n * 


(A-3) 


<A-4) 


APPENDIX  B 


DERIVATION  CP  THE  da/dN  EQUATION  PGR  THE  QUADRATIC  LOG-LOG  T'POINT 
INCRB»SNTAL  FOLYNOHIAL  WTUOD 

Th«  fitted  polynoAlal  equatioa  for  the  quadratic  log- log  7-polnt  In- 
creaental  polynomial  mathod  la  given  by 


logic  • • **o  ^ \S  ^ ’’2  "lS 

where  Ni^  la  given  by  equation  A-2.  Substituting  into  Equation  B-1  for 


' logic  ” ‘ 

^0*10  • • *>o  j ^ ^2 


Letting 


c . l!!iiLiLlfi 


(B-2) 


(B-3) 


Then 


Solving  for  a, 


^ L 


dK  N . In(lO)  . C- 


‘»*10  * • ".  * ‘i  " * ‘2 


b + b,  U + b- 
o 1 2 


a - 10 

Taking  the  derivative  of  a with  respect  to  U, 


(B-A) 


(B-5) 


(B-6) 


Ifl 


10 


• 10 


10 


u 

du 


J - 1^10  ‘ . 10  ■ 2 bj  U In  (10) 


b,  b U . 

+ 10  ^ • 10  ^ . bj  ln(10)J 


(B-8) 


Then 


4§  - 10  -10 


b,  U bj 

• 10 


In  (10)  • , 2 b.  U + b,  : (B-9) 

u 2 Ij  . . , 


Ceix:  the  cheln  rule, 


. 4a  . is 

dN  dU  dN 


(B-10) 


bo  bj  U b^U^ 


- 10  *10 


10 


• ««  • L“j  “ * V • L«~T..aor:ig. 


1>„  b,  U b.O* 
' o 1 Z 


10  • 10 


10 


<2b^  U b^) 


(B-ll) 


Cj  • H 


Subetlcutlng  for  U end  evelueelng  et  the  Bidpelnt,  R^, 

b.lo»R.  b.C.^  b,  (logR,)*  - ZbaCjlogMj  + bjC^^- 


APPENDIX  C 


KRIVATIOM  OP 


From  linear  regression,  the  coefficient  of  multiple  determination, 


R^,  Is  given  by  [34] 


R‘  - 1 - 


SSRBS 


(C-1) 


Where  SSRBS  Is  the  residual  sun  of  squares  and  TCSS  Is  the  total  corrected 
sum  of  squares.  Since  SSRBS  and  TCSS  are  srasured  In  the  vertical  dlrec* 
tlon,  It  was  desirable  to  correct  then  so  that  their  direction  is  normal 


to  the  slope,  m.  Let  the  slope  be  given  by 


(C-2) 


Where  J Is  the  side  of  a triangle  along  the  least  squares  line  and  K is 
Che  side  of  the  triangle  perpendicular  to  the  least  squares  line.  From 


basic  geomstry. 


(C-3) 


Where  1 la  Che  third  side  of  the  triangle.  I Is  always  in  a vertical  dl* 
rectlon.  Substituting  from  equation  C-2  into  C-3  for  J, 

- (n  . K)^  + 


2 2 

- K (.n  + 1) 


(C-4) 


Solving  for  K , 


m*  + I 


Solving  for  the  correction  of  the  elope, 


(C-5) 


7 


(C-6) 


+ 1 


Plugging  chle  into  equetion  C-1  to  obteln  k corrected  for  the  elope. 


called  C , 


C - R 


(K^/1^) 

SSRES~,  1 ' 

^ + i- 


(C-7) 

(C-8) 
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APFBMDIX  D 


DNDDPG  DOCOMRiUAIIOM 

Thxs  program  conalscs  of  a main  program  and  25  subreuclnaa.  Tha  main 
program  (DNDDPG)  reads  in  the  daairad  ^a  vs.  AM  data  and  calls  subroutine 
KLTA  to  rs'create  the  original  a va.  N data.  The  program  flow  is  then 
transferred  to  subroutine  CLASS  which  divides  the  data  into  constant  Aa 
data  sets  and  then  calculates  the  histogram  frequencies  for  each  constant 
Aa  data  set.  The  program  flow  is  then  transferred  to  subroutine  STPLOI 
which  determines  the  distribution. 

Subroutine  STPLOT  uses,  directly  or  indirectly,  the  following  sub* 
routines. 

1.  Parameter  Bstimation  Subroutines 

1.  OOIi«N 

2.  cRvprr 

II.  Scaling  Subroutines 


1. 

FRBFLr 

2. 

WBLF&B 

3. 

LGSCAL 

4. 

LN8CAL 

3. 

iNuec 

6. 

ODSCAL 

III.  Plotting  Subroutines 
1.  DOFLOT 
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2. 


MRMPLI 


3.  U)6FLT 

4.  miPLI 

5.  OOiaiS  -.i 

IV.  Output  Subrout Liws 

1.  RIIMX 

2.  RITPAR 

y.  fioQoral  PurpoM  SubrouciMB 

1.  LSTSQR 

2.  BANK 

3.  OUTLIR 

4.  »MIAB 

5.  SIMPSN 

6.  MAXR 

7.  MAXI 

A listing  of  this  program  can  be  obtained  from; 

Prof.  B.  M.  Hillberry 

School  of  Mechanical  Engineering 

Purdue  University 

West  Lafayette,  Indiana  47907 

Phone  (317)  494-1600 


APPENDIX  K 


CCDDP  DOCOWNTATIOM 

This  program  consists  of  a main  program,  53  subroutines,  and  18  func- 
tion subprograms.  The  main  program  (CCDDP)  reads  in  the  desired  replicate 
cycle  count  data  and  writes  it  by  calling  subroutine  &TTDAT.  The  program 
flow  is  then  transferred  to  subroutine  CLASS  which  calculates  the  histo- 
gram frequencies  for  the  data.  The  program  flow  is  then  passed  to  sub- 
routine STPLOI  which  determines  the  distribution. 

Subroutine  STFLOT  uses,  directly  or  indirectly,  the  following  sub- 
routines and  function  subprograms. 

I.  Parameter  Estimation  Routines 

A.  Subroutines 

1.  MLBLN 

2.  MLEW 

3.  MLEG 

4.  KLEGG 

5.  HJ 

B.  Supporting  Function  Subprograms 

1.  PLN 

2.  FW 

3.  FG 

4.  FOG 
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1 
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j 

i 


i 


i 

1 

i 


i 
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n.  Statistical  Paranatara  Subroutinaa 

1.  WtSTAT 

2.  WBSTAT 

3.  GMSTAT 

in,  Goodnaaa  of  Pit  Routinas 

A.  Subroutinaa 

1.  CHISQR 

2.  KOLS)« 

3.  MRMCS 

4.  HBLCS 

5.  GAMCS 

B.  Supporting  Function  Subprograms 

1.  PNRH 

2.  FWBL 

3.  FGAM 

IV.  Output  Subroutines 

1.  RITPAR 

2.  RITRBS 

3.  PAROUT 

V.  Plotting  Routinas 

A.  Main  Plotting  Subroutinaa 

1.  AVNPLT 

2.  HISPLT 

3.  NRMPLT 

4.  LOCPLT 

5.  WBLPLl' 

6.  GAHPU 
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B.  Supporting  Plotting  Subroutinea 

1.  ODAXIS 

2.  LGAXIS 

3.  GKAXZS 

4.  SCIMOr 

VI,  Scaling  Subroutinea 

1.  lAMSCL 

2.  WBLSCL 

3.  GAMSCL 

! 4 . LGSCAL 

5.  LMSCAL 

i 6.  INLNSC 

7.  ODSCAL 

8.  SCAUL 

VII.  Streaa  Intensity  Calculation  Routines 

) 

A.  Subroutine 
1.  DELTAK 

i 

B.  Function  Subprogram 
1.  FAB 

Vlll.  General  Purpose  Statistical  Routines 

A.  Subroutines 
1.  NRMTAB 

^ 2.  omiR 

B.  Function  Subprograms 

I 1.  PNORM 

^ 2.  FGAMIA 
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4 


FTRIGK 


5.  FIKGAM 

6.  FGM 

7.  FGMmr 

8.  FGMNBG 

9.  FSER 

10.  FFRAC 

IX.  General  Purpose  Subroutines 


1. 

TABL£L 

2. 

INTUAV 

3. 

IMTBRP 

4, 

INVINT 

5. 

LSTSQR 

6. 

INVMAT 

7. 

RANK 

8. 

HANCHA 

9. 

INITR 

10. 

ITOR 

11. 

LOG 

12. 

HAXR 

13. 

HAXl 

A listing  of  this  program  can  be  obtained  from  Professor  B.  M. 
Hlllberry  (Appendix  D). 


APHIIDIX  P 


C681H>P  DOCtnCNTATION 

Tbit  progrta  coatltet  of  • ntln  progrMB,  SO  aubroutloet,  and  17  func- 
tion aubprograat.  This  prograa  is  naarly  Identical  to  the  CCDDP  prograa 
(Appendix  B)  and  only  the  main  program  (CORSDP)  and  3 tubroutines  are 
changed.  These  3 aubroutioes  are; 

1.  CLASS, 

2.  STPLOT,  end 

3.  RITDAT. 

This  prograa  requires  an  input  of  replicate  growth  rate  data  and  has  the 
saae  output  as  the  CCDDP  program.  Subroutines  DELTAK,  ITOR,  and  HAXI  and 
function  subprogram  FAB  need  not  be  loaded  for  this  program.  A Hating 
of  this  program  can  be  obtained  from  Professor  B,  M,  Hillberry  (^pendix 
D).  ^ . 
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AFRMDIX  G 

DION)?  DOCUICNIAfZON 

This  program  coosliCg  of  a main  progrM,  49  subroutlnaa,  and  17 
function  aubprograma.  This  program  la  naarly  Idantical  to  tha  C(Slfl>P 
program  (Appandix  F)  and  only  tha  main  program  (DHDDF),  9 aubroutinaa, 
and  3 function  aubprograma  ara  changed*  Tha  9 aubrouClnaa  that  ara  >- 
changed  ara ; 

1.  CLASS, 

2.  8TPL0T, 

3.  MUG, 

4.  MUGG, 

5.  GWTAT. 

6.  GANGS, 

7.  SZTBlAT, 

8.  RITFAR,  and 

9.  OAMPLT. 

Tha  3 function  aubprograma  chat  ara  chaagad  ara; 

1.  FC, 

2.  FOG,  and 

3.  FOAM. 

Thia  program  raquiraa  an  Input  of  rapllcata  growth  rata  data  and  haa  tha 
aamn  output  aa  tha  CGBDDP  program,  lubrouclna  IHVNAI  noad  not  ba  lomdad 
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tor  thi.  ,ro,t...  A llottn,  of  thl.  ,ro,r«.  b.  obtolood  fr« 
Profeoor  B.  M.  Hillb«rry  (Appandix  D). 
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AFFBHDIX  H 


OSLTCF  DOCiaeiTTAIIOIf 

This  progrsm  consists  of  s miiln  program  and  5 subroutlnas.  Tha  main 
program  (DELTCP)  raads  In  the  desired  a vs.  N data  and  calls  the  proper 
subroutine (a)  to  calculate  the  Ae  vs.  AN  data  according  to  the  desired 
calculation  method  chosen.  The  Aa  vs.  AN  calculation  subroutlnas  are; 

1.  REMOVE. 

2.  STRIP,  and 

3 . DELTA. 

The  main  program  then  calls  subroutine  RITDAI  to  write  the  A*  vs.  AN 
data.  The  only  general  purpose  subroutine  required  Is  subroutine  LST8QR. 
A listing  of  this  program  can  be  obtained  from  Professor  B.  M.  Hlllberry 
(Appendix  D). 


2M 


APRNDZX  X 


OAONCP  OOCiaSNIATKM 

Thl«  program  coaalata  of  a main  program,  22  aubroutlnaa,  and  1 func- 
tion aubprogram.  Tha  main  program  (DADNCP)  raadi  in  tha  daairad  a va.  M 
data  aet  and  calls,  diractly  or  indiractly,  tha  following  subroutinas  and 
function  subprogram, 

I.  Growth  Rata  Calculation  Subroutinas 

1.  DADN 

2.  SECANT  

3 . M0D3EC  '■ 

4.  STRIP 

5.  BVAL 

IX.  Error  Oatarmlnation  Subroutinas 

1 . DELTA 

2.  IIITECR 

3.  ERROR 

III.  Output  Subroutinas 

1.  RITDAT 

2.  RITSBS 

3.  RISVIT 

XV.  Plotting  Subroutinas 
1.  AVMPLT 


2 


LDPLT 


V,  ScAlias  St^routiiw 


1.  USCAL  , 

Stress  Intaniltf  Calculatlea  liMtlM* 

A.  SnbwtlM 
1.  VD/m 

B,  Fuoetloa  8u^ros<«a 


1.  lAB 


. VII.  OcMMl  Purpose  Subroueinee 


£?rl.  inm 

2.  IKITI 

3.  CHBCK 

4.  ITCa 

5.  MAXI 

6.  LOO 

7.  LST8QR 


A Ileclng  of  thle  progren  can  be  obtained  from  Profeaaor  B.  M.  Hlllberry 
(Appendix  D) . 
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AVNFRD  DOCtnCHTATlQN 

This  program  consists  of  a main  program,  19  subroutines,  and  8 func- 
tion subprograms.  The  main  program  (AVHnB)  reads  in  the  desired  distri- 
bution parameters  and  calls,  direetlp  or  indirect^,  the  following  sub- 
routines and  function  subprograms. 

1.  Random  Number  Generating  Subroutine 

1.  RMGEM 

11.  Inverse  Distribution  Subroutines 

1.  IHVD18 

2.  INVNRM 

3.  1MV21M 

4.  1HV3LH 

3.  INWBL 

6.  INVCAM 

III.  Prediction  Subroutines 


V.  Plotting  Subroutlnn 


1.  AVNPLT 

VI.  Gtnaral  Purpos*  Statistical  Butinas 

A.  Subroutine 

1.  ISMTAB 

B.  Function  Subprograms 

V ..  ^ PNORM 

2.  FGAMtA 

3.  PtNGAM 

4.  FGM  - : 1.  - 

5.  FCMINT 

6.  FGMNBG 

7.  FSBR 

8.  FFRAC 

VII.  General  Purpose  Subroutines 

1.  TABLEL 

2.  IKTHAV 

3.  IinBRP 

4.  iNvnrr 

5.  MAXR 

A listing  of  this  program  can  be  obtained  from  Professor  B. 
(Appendix  O). 
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APRHDIX  K 
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The  rando«iM4 
fftlJUsw*. 


XAMpON  non  OF  EXRRDSMEAL  TSSTS 

•rd*r  oi  tb|  68  ApcciMBS  ua«d  during  testing  is  ss 


Vt  .v.  y- 


■ . « ..  - 

- 1. 

U1 

49.  ‘ 

141 

- 37. 

49 

55. 

127 

101 

20. 

60 

38. 

111 

56. 

6 

1*7 

—^21. 

85 

• ^ ^39. 

102 

57. 

71 

•-  -4. 

77 

22. 

50 

40. 

64 

58. 

58 

66 

134 

-82 

59. 

78 

■■■  ■ 

- 

t tty 

•• 

. 

-'"i. 

119 

24. 

55 

42. 

94 

60. 

106 

■ ■ ^,7. 

*7 

23. 

14 

•743. 

138 

61. 

129 

112 

-26. 

21 

44. 

130 

62. 

31 

»ygw<9 

- ...  ‘r^t 

. . T ‘ 

■ - ■ -^  -.x-9:- 

*2 

- .-27. 

7 

45. 

143 

63. 

125 

10. 

40 

28. 

2 

46. 

57 

64. 

32 

11. 

18 

29. 

70 

47. 

41 

65. 

116 

12. 

132 

30. 

19 

48. 

9 

66. 

62 

13. 

118 

31. 

32 

49. 

74 

67. 

14. 

1 

32. 

83 

30. 

33 

68. 

121 

13. 

35 

33. 

149 

51. 

144 

16. 

123 

34. 

80 

52. 

81 

17. 

37 

35. 

97 

53. 

93 

18. 

139 

36. 

96 

34. 

99 
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